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EXECUTIVE SUMMARY 


An experimental and theoretical study of radiative feedback, 
burning rates, radiative heat loss fractions, and flame heights 
for pool fires with diameters ranging from 4.6 cm to 100 cm was 
completed. Transient measurements of soot and temperature 
distributions were obtained in 7.1 cm and 30 cm heptane and 
toluene fires using a three-wavelength emission/absorption probe. 
The heat release rates of the fires varied from 0.6 kW to 2166 kW 
allowing a study over a wide range. A variety of fuels were tried 
but most of the measurements were restricted to methanol, heptane, 
and toluene as representatives of the alcohols, paraffins and 
aromatics. Radiative feedback was measured using a new insitu 
purged optical probe inserted at the level of the fuel surface. 
Measurements of reflection of energy from the fuel surface were 
also obtained. A multi-ray radiation calculation procedure 
utilizing simultaneous single-point measurements of soot volume 
fractions and temperatures was used to estimate the heat transfer 
to the surroundings and that to the fuel surface. Importance of 
turbulent fluctuations of different frequencies on the radiation 
heat flux was studied using filtered simulations. Important 
conclusions of the present study are: 


(1) The total radiative heat loss fractions of several alcohol 
fuels were approximately 0.2 and those of paraffin and 
aromatic fuels were approximately 0.3 within the range of 
pool diameters and heat release rates, 


(2) Flame heights of alcohol and paraffin fuels were in 
reasonable agreement with existing correlations while those 
of aromatic fuels showed possible effects of lower combustion 
efficiency, 


(3) Burning rate measurements showed significant effects of fuel 
type for all pool diameters. For aromatic fuels such as 
toluene, the radiation dominated limit was reached at a 
relatively small pool diameter (between 5-7 cm). For heptane 
fires, the limit was reached at approximately four times the 
value for toluene. However, none of the pools in the present 
size range showed a clear conduction regime. The experiments 
at NIST with much smaller diameter pool flames conducted by 
Dr. John Yang should shed further light on this. 


(4) Radiative heat flux incident on the fuel surface was large 
(almost equal to the total heat feedback) for all the fuels 
for pool sizes above 30 cm near the pool center. For toluene 
fires, the heat flux remained high and almost constant as a 
function of radial distance. For heptane fires, the 
radiative feedback diminished a little with radius. However, 
for methanol fires, a significant decrease was observed with 
radius. It was found that the burning rate vs. pool diameter 
drop showing regimes of different heat transfer mechanisms is 
not consistent with insitu radiative heat feedback 
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(6) 


(7) 


measurements. Hence, the utility of such a map for fire- 
modelling should be reconsidered. 


Measurements of transient emission and absorption properties 
resolved to 10 mm show that: 


e Large fluctuations in emission intensity and soot volume 
fractions occur in all regions of the flames. 


* A layer of relatively cold soot particles exists near the 
fuel surface probably resulting from hydrodynamic transport 
from the edge of the fire. 


e Comparison of measurements of soot volume fractions based 
On absorption with those based on emission shows that 
within a 10 mm long probe volume, a large portion of the 
soot particles are at relatively low temperatures and do 
not contribute to the emission. 


Calculations of radiation intensity and radiative heat flux 
using the time series of emission soot volume fractions and 
temperatures showed that: 


¢ The 10 mm resolution offered by the present probe was 
sufficient for obtaining effective values of these 
properties. 


¢ Use of mean values of temperature and soot volume fractions 
May result in an order of magnitude underestimate of 
radiation intensity and heat flux. However, due to the 
high sensitivity to temperatures, the error is normally 
corrected by using slightly higher average temperatures. 


¢ Use of a Gaussian filter applied to the temperature and 
soot volume fraction time-series showed that relatively low 
frequency fluctuations in temperature and soot volume 
fractions contribute the most to the increase in radiation 
heat flux above that estimated from mean properties. The 
frequency above which the fluctuations appear to have 
minimal effects is close to the characteristic oscillation 
frequency of the fires. 


Predictions of radiative heat feedback and heat flux to 
surroundings using a multiray method showed that: 


¢ The total heat flux changed by less than 15% between 22 and 
450 rays. 


e Reasonable comparison between predictions and measurements 
was observed. 


The multi-ray procedure allows an examination of the regions 
of origin of the radiative heat feedback. It was observed 
that bulk of the feedback originated at heights below 0.5 D 


in toluene fires, and at heights below 3.5 D in heptane 
fires. 


The results of the present research are being applied by NIST 
personnel to the development of advanced global burning rate 
models. The multi-wavelength emission absorption probe has been 
used and improved by NIST researchers. The measurements have 
spurred activity in the areas of optical properties of soot, flame 
structure that leads to the relatively low temperature soot, 
utility of absorption vs. emission measurements, use of 
appropriate effective temperature from the fire, and evaluation of 
current entrainment rate correlations and data for compatibility 
with global fire models. Improvements in such models resulting 
from the knowledge gained in studies, such as the present work, 
are expected to contribute to long-term fire safety. 
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Chapter I 
Introduction 


1.1 General Statement of the Problem 

Radiation from liquid fueled flames is an important heat 
transfer mechanism in combustion systems such as fires, 
furnaces and gas turbine combustors. Radiative heat flux from 
fires to the surrounding objects determines the possibility of 
ignition, flame spread and flashover. Burning rates of 
objects on fire determine safe egress times, heating rates of 
surrounding objects and flame spread rates. Knowledge of the 
mechanisms that determine the burning rate and radiative heat 
flux of fires has been sought since early days of fire science 
(Hottel, 1958). 

GeRis (1979), Mudan (1984) and Quintiere (1992) have 
reviewed past work concernging burning rates and radiation 
properties of pool fires. Existing methods of calculating 
radiation from liquid fueled fires are based on homogeneous, 
isothermal flame approximations and an effective average flame 
temperature and emissivities. The effect of the differences 
in the structure of the flames generated by different fuels 
are considered only in a global manner in this approach. 
Effects of turbulent fluctuations on radiative transfer are 
not explicitly considered. Motivated by these observations, 
the overall objective of the present investigation was to 
study the radiation properties of liquid fueled diffusion 
flames. 

The study concentrated on axisymmetric diffusion flames 
from circular burners containing liquid fuels. The main 
research issues addressed here are: (1) total radiative output 
of liquid fueled flames; (2) radiative heat transfer to the 
fuel surface; (3) local emission and absorption properties of 
strongly radiating pool flames; and (4) the significance 
turbulence-radiation interactions in buoyant’ turbulent 
luminous pool flames. Particular emphasis was placed on the 
mechanism of radiative heat transfer to the fuel surface, 
termed radiative heat feedback, in luminous flames where the 
continuum radiation properties depend on local soot volume 
fractions and temperatures. 

The burning rate of pool fires is determined by the rate 
at which energy is transferred to the liquid fuel. Fig. 1-1 
shows the energy balance at the fuel surface of a liquid pool 
fire. PH is the energy from the flame incident on the fuel 
surface, od +5 the energy conducted into the fuel by the 
burner varies mH,, is the total latent heat of vaporization, 
and Q.ind-loss 1S the energy conducted through the fuel layer. 
The energy feedback term Q, involves radiation, convection, 
and conduction from the hot flame gases to the liquid fuel 
surface. This term proves to be difficult to characterize 
since turbulence and radiation are often involved in the 
mechanisms. 

As discussed in the literature reviews cited above, the 
related studies cited in the next section and the literature 
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reviewed in the introductory sections of each of the following 
chapters, there are specific breachs in our understanding of 
the effects of fuel properties, burner size and sooting 
tendency on heat transfer from the flames to the liquid 
surface (heat feedback). Information concerning the relative 
contribution of radiative and convective processes to the 
heat feedback is’ lacking. Furthermore, the effect of 
transient local flame structure (distribution of scalar 
properties such as species concentrations that affect 
emissivity and temperature) on heat feedback has not be 
quantified. In situ measurements of the heat feedback and the 
local, transient scalar properties are lacking for many fuels 
of interest. 

The present work addresses some of these gaps. The main 
technical application of this research is to develop a 
fundamental understanding of the continuum’ radiation 
properties in luminous diffusion liquid fueled flames, 
particularly with regard to the mechanism of heat feedback. 
The results have application to modelling fires within 
structures, material characterization and test methods and 
fire detection systems. 


1.2 Related Studies 

Past work in radiation from pool fires has been reviewed 
by deRis (1979), MudQn (1984) and Fischer (1988). These 
reviews provide descriptions of the approximate radiation 
calculations methodologies and experimental data concerning 
pool fires. Fernandez-Pello (1991) and Quintiere (1992) have 
reviewed the fundamental mass and energy transfer processes in 
pool burning. Faeth et al. (1989) and Fischer (1988) have 
reviewed the effects of fluctuating scalar properties in flame 
radiation. Background material on radiation in pure and 
particle laden gases is provided by Siegal and Howell (1981), 
Ludwig et al. (1973), Goody (1964) and Hottel and Sarofim 
(1967). 

Burgess and Hertzburg (1974) and Blinov and Khudiakov 
(1957) have shown that for certain fuels, the radiative 
feedback is dominated by conduction below a pool size of 10 cm 
diameter, between a pool size of 30-50 cm diameter convection 
is important and above this range the fires become radiatively 
dominated. For large fires, the relative contribution of 
radiative and convective heat flux may also depend on the 
position on the pool surface due to changes in flame shapes 
and convective velocities. 

In non-luminous flames, energy is radiated from gases 
such as CO,, CO and H,0. Edwards and Menard (1964) developed 
wide band exponential models to predict the emission from 
radiating gases which has been used successfully by Negrelli 
et al. (1977) and Liu et al. (1981) in laminar methane/air 
flames. Grosshandler and Sawyer (1978) examined the non- 
luminous spectral radiation properties of methanol/air 
combustion products. This study used a statistical narrow- 
band model to make encouraging predictions in the post-flame 
region. 


Several recent studies have made transient measurements 
of radiation and flow field properties of pool fires. The 
results of these studies indicated the fluctuating nature of 
the scalars in all regions of pool fires. McCaffrey (1979) 
measured velocity and temperature distributions in natural gas 
pool fires. This study identified the various flow regions 
within a pool fire. Weckman and Sobiesiak (1989) studied the 
entrainment mechanism in medium scale acetone pool fires. 
Bouhafid et al. (1989) measured temperature, gas species 
concentrations and monochromatic absorption coefficients near 
the fuel surface of kerosene pool fires. They identified the 
fluctuating nature of scalars in pool fires. Fischer and 
Grosshandler (1987) measured time-averaged compositions of 
radiating gases and soot, and fluctuating temperatures in 
alcohol pool fires. Shinotake et al. (1985) made transient 
measurements of heat feedback in heptane pool fires. 

The effects of scalar fluctuations on luminous radiation 
have been shown to be significant. Mean property methods that 
are frequently used to compute radiation intensities 
underestimate the mean radiation intensity. This has been 
noted by Cox (1977), Faeth et al. (1988), Gore and Faeth 
(1986, 1987), Kabshinov and Kmit (1979) and Portscht (1974). 
Stochastic simulations following the time series analysis 
methods of Box and Jenkins (1976) have been successfully 
applied to fully turbulent jet flames by Kounalakis et al 
(1988a,b). Sivathanu et al. (1991) and Sivathanu and Gore 
(1991) have shown that the joint statistics of temperature and 
soot volume fractions are necessary for obtaining reasonable 
predication of radiation intensity in strongly radiating 
flames. 


1.3 Overall Objectives 

The present study seeks to improve the understanding of 
radiation from liquid fueled pool flames by addressing some 
open research issues. These include: 


1. The effect of fuel type and scale on mass burning 
rates and flame radiation; 

2° The mechanism of radiative heat transfer in pool 
fires; 


3. Local emission/absorption properties of luminous pool 
fires; and 

4.  Turbulence-radiation interactions in prediction of 
radiative feedback in luminous pool fires. 


Chapter II describes the global properties of all the 
flames studied in the present investigation. The global 
properties include mass burning rates, flame heights and 
radiative heat loss fractions. Chapter III describes the 
measurements of radiative heat feedback to the fuel surface in 
medium sized pool flames. Chapter IV reports the measurement 
of local emission/absorption properties in a small scale 
strongly radiating pool flame. The fluctuating nature of 
emission scalars is investigated. Chapter V describes 


measurements of emission/absorption properties in a medium 
scale pool fire and the use of these properties in the 
prediction of radiative heat feedback. 


CHAPTER II 
GLOBAL FLAME PROPERTIES 


2.1 Introduction 

Pool fires are generally luminous and emit radiation due 
to both soot and combustion products. The energy radiated 
from the flames is an important parameter in fire safety. The 
amount of heat flux directly contributes to the hazard posed 
by the flame and the spread rate of the fire. 

A number of experimental studies have considered the 
total radiant heat loss from gaseous turbulent and laminar 
diffusion flames (Markstein, 1984; deRis, 1979; Faeth et al, 
1989; Gore, 1986; Sivathanu, 1990; and Gore et al, 1991). 
Markstein (1984) found that the radiative heat loss fraction 
is a constant for highly buoyant turbulent flames, as opposed 
to laminar gaseous flames for which the radiative fraction 
increases until the smoke point (defined as the fuel flow rate 
at which visible dark smoke is just emitted at the flame tip) 
is reached. Jeng and Faeth (1984) report the radiant heat 
loss from methane flames both in buoyant and in the forced 
regimes. Subsequent studies (Gore, 1986; Sivathanu, 1990) 
support the experimental finding that for low speed turbulent 
flames burning a variety of fuels (acetylene, carbon monoxide, 
methane, ethylene, propane, propylene, and ethylene), the 
radiation losses are a constant fraction of the chemical 
energy release. 

Past measurements of radiant heat loss from liquid and 
solid fueled flames have generally relied on single location 
measurements. This method assumes that the radiative heat 
flux from flames is spherically isotropic: 


O-aq7 4% R*qQ" (2-1) 


where R is the radius of the sphere and q" is the radiative 
flux at the measurement location on the sphere. Modak and 
Croce (1977) used this method for estimating the radiative 
heat loss fraction from PMMA. Koseki and Yumoto (1988) 
measured the radiant heat loss for large scale heptane fires 
relying on single point measurements. Bouhafid et al. (1989) 
measured radiative heat flux at various distances from pool 
fires and calculated the heat loss fraction. Recently, Kédylu 
and Faeth (1991) have reported measurements of radiative loss 
fractions for 0.1 (m) diameter flames stabilized on a porous 
matrix. 

The objective of this portion of the investigation was to 
determine the total radiative output of liquid-fuel pool 
flames at a range of scales burning a variety of fuels. Faeth 
and co-workers (Jeng and Faeth, 1984; Gore, 1986; and 
Sivathanu 1990) utilized a multi-location method for radiant 
heat flux measurements in gaseous diffusion flames. This 
technique, used in the present investigation, avoids the 
assumptions used in single point measurements, and had not 
been utilized in previous studies of liquid and solid fueled 


flames. 

Additionally, measurements of effective temperature and 
transmittances were made for two flames of different sooting 
tendency. Measurements of the spectral flame emission and 
spectral fuel vapor transmittance were made for heptane and 
toluene. Measurements of flame heights and fuel burning rates 


were also made. 


2.2 Experimental Methods 
2.2.1 Experimental Apparatus 

The experiment involved the burning of liquid fuels from 
a circular pan in a quiescent room environment. The flames 
were operated in a large enclosure, the location and size of 
which varied due to the large range of burner diameters which 
were used in the study. A sketch of the - general 
characteristics of the enclosures is shown in Fig. 2-1. The 
enclosure utilized wire mesh screens and/or plastic sheets to 
minimize effects of room disturbances on flame shape and other 
flame properties. Each facility employed an adjustable 
exhaust system to remove combustion gases and soot from the 
enclosure. The burners were mounted rigidly and was placed 
approximately 15 cm above ground level. 

Pool burners of inner diameter 4.6, 7.1, 30, 60, 100 cm 
were used. Other relevant dimensions of the burners are given 
Table 2-1. The 4.6 and 7.1 cm burners (Fig. 2-2) were 
uncooled and constructed from pyrex and included ports from 
which fuel could be delivered and a thermocouple could be 
installed. 

The larger burners were made of rolled steel and are 
depicted in Fig. 2-3. Each burner contained several ports in 
the base for delivery of fuel and installation of thermocouple 
arrays. A cooling water reservoir was placed on the bottom of 
the steel pan and was fed cold tap water throughout each test. 
Incoming and out-going cooling water temperatures were 
monitored throughout the test by chromel/alumel thermocouples 
connected to a data acquisition computer. The thermocouple 
arrays consisted of individual thermocouples which were placed 
at different distances relative to the burner lip. The output 
of these thermocouples provided a temperature profile of the 
fuel in the liquid phase. 

Fuel burning rate was monitored using a load cell mounted 
under the fuel supply reservoir. In a few trials with the 
smaller pools (4.6 and 7.1 cm), a 75 ml burette and stop watch 
were used and the burning rates manually recorded. The output 
of the load cell was collected by a data acquisition computer. 
The uncertainty of the load cell data is estimated to be 5% 
based on repeated measurements and calibrations. The 
uncertainty of the burette data is estimated to be 10%. 

Preparation for the tests began with activation of the 
cooling water. Sufficient time was allowed for the cooling 
water temperature to stabilize. The burner was filled until 
the fuel level was approximately 1 cm below the edge. The 
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fuel reservoir was filled and any air in the fuel line was 
bleeded out. The regulating and solenoid valves were shut and 
exhaust system activated. After ignition, the exhaust system 
was adjusted so that combustion products were removed but the 
flame shape was unaffected. The fuel regulating valves were 
adjusted to maintain the fuel at the required level. 

The fuel level was maintained constant relative to the 
burner lip during the entire test. The level was 5 mm below 
the edge for the larger burners (D = 30 cm) and 3.5 mm below 
the edge for the smaller burners. The level control was 
accomplished using a gravity feed fuel system which utilized 
a fine adjustment regulating valve and a temperature- 
controlled solenoid valve. The control system is depicted in 
Fig. 2-4. The solenoid valve controller received input from 
a thermocouple placed just above the liquid fuel surface. The 
solenoid valve is activated open when the temperature rises 
above a prescribed value and shuts when the temperature limit 
is again reached. The regulating valve controls the rate at 
which fuel enters the burner, preventing (along with the 
baffles placed at the bottom of the burner) disturbances of 
the fuel surface by the incoming fuel. 

Data collection began only after a warm-up period of at 
least 10 minutes. This allowed for the fuel in the pan and 
the surroundings to attain a relatively steady temperature. 
Obtaining acceptable flame shapes (negligible amount of flame 
lean) became more difficult as the scale increased. This was 
true even with the enclosures described previously. Much 
effort was made to attain axisymmetric, uniformly attached, 
vertical flame shape with minimum room disturbances during the 
warm-up period. The effort involved placement of screens and 
baffles and adjustment of exhaust flow rate to achieve 
desirable conditions. In the larger scale tests (D > 30 cm), 
the fuel height in the burner was varied slightly (5 - 15 mm) 
in order to maintain reasonable flame shapes. This was 
necessary for heptane and toluene as fuel vapor would 
occasionally ignite in the over flow ring on the burner if the 
fuel level was at 5 mn. 

2.2.2 Instrumentation 

Flame heights were obtained from video tapes (30 
frames/sec) and still photographs. The luminous flame heights 
were obtained by averaging approximately one hundred frames 
selected at random. When necessary and possible, room 
lighting was dimmed and dark backgrounds were used to enhance 
the quality of the video and stills. Distances on the film 
and video where calibrated by photographing’ scales. 
Experimental uncertainties for these measurements depended on 
the room lighting, flame luminosity and film sensitivity. For 
a specific set of these conditions, the uncertainties were 
less than 10% but could be as high as 40% for different tests. 

A calibrated radiative heat flux transducer was used to 
measure total radiative heat fluxes. The method used to 
determine total radiative heat flux is identical to that Gore 
(1986) and Sivathanu (1990). The transducer(s) used were 
water cooled and provided a 120° to 150° viewing angle. The 
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Fig. 2-4 Sketch of the fuel supply and regulation system 
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radiant heat flux incident on locations on the burner plane 
and also on locations parallel to the flame axis was measured. 

In the burners with diameters less than 30 cm, a single 
heat flux transducer was traversed both radially and axially 
to capture the radiant energy. In the larger sized pools, 
test times were generally much shorter due the large amount of 
fuel consumed. In these tests, an array of transducers was 
used. The transducers were fixed in place in both the radial 
and axial directions. The transducers were calibrated using 
a known standard transducer and a tungsten lamp source. 
Experimental uncertainties for the total radiant output 
measurements are estimated to be less than 15 percent based on 
repeated tests. 

Effective temperatures and transmittances were measured 
by using three wavelength optical pyrometry. A schematic of 
the experimental apparatus is shown in Fig. 2-5. The 
technique consists of the simultaneous measurement of flame 
emission intensities at two wavelengths (A = 0.9 and 1.0 mum) 
and laser transmission (A = 0.632 um) across the entire width 
of the flame. The light leaving the flames is collected by a 
0.6 cm diameter, 25 cm long stainless steel tube that has a 
nitrogen purge of 1 cc/min. A second tube with identical 
purge flow serves as the cold background for the measurement 
and also as the guide for the incident He-Ne laser used for 
the transmittance measurements. Two beam splitters are used 
to divide the light into three parts. Focussing lenses and 
line filters with 10 nm bandwidth on each channel complete the 
optics. The emission detectors are two photomultipliers tubes 
and the transmittance detector is a laser power meter. The 
incident laser power is monitored by an identical detector. 
All signals are sampled by a laboratory computer at 250 Hz. 

Spectral fuel vapor transmittance is measured using the 
experimental apparatus shown in Fig. 2-6. Two probes, 
defining a 10 mm path length, are placed at a position in the 
flame were little emission and soot is present (x/D = 0.01). 
Blackbody radiation tranmission is monitored through the 
probes. A 1/8 m grating monochromator is used to determine 
the spectral content of the transmission signal. A Pb-S or 
Pb-Se detector is utilized depending on the wavelength of 


interest. The detection signal is sampled by a laboratory 
computer. 

Spectral flame emission is measured for a path length 
that encompasses the entire width of the flame. A 1/8 m 


grating monochromator defines the wavelength interval of 
interest and a Pb-S or Pb-Se detector is used for the 
measurement. An apporpriate optical edge filter is used to 
prevent measurement of harmonics. The signal is monitored by 
a laboratory computer at 250 Hz. The detectors are calibrated 
at each wavelength by a blackbody radiation source. 

Fuel surface reflection of incident infrared radiation 
(A= 1.0 - 3.0 um) is measured using the apparatus depicted in 
Fig. 2-7. An alignment laser is used to position the mirrors 
and infrared detector(s) and measure the incident radiation 
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Fig. 2-7 Apparatus for measurement of fuel surface reflectivity 
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angle. Infrared radiation from a blackbody source is 
colimated used 6 mm diameter tubes. The detector signal is 
measured with the use of a lock-in amplifier and a laboratory 
computer. 

2.223 Data Analysis 

The radiative heat flux was measured in a direction 
normal to the plane of the fuel surface at the burner plane 
for a distance R and then in a direction perpendicular to the 
flame axis (see Fig. 2-4). These measurements where extended 
in the axial direction until close to background flux were 
found. Radiation incident on the pool surface were assumed 
equal to value measured just outside the burner. The radiant 
flux distributions are integrated to enclose the flame in a 
"Semi-infinite" cylinder (see Hamins et al. (1991)), thus 
providing the total radiative energy (Q, aq) of the flame. The 
radiative heat loss fraction (X,) is defined as: 


Oy ag hom (2=2) 


where Q,., is the fuel burning rate (m,) multiplied by the 
ideal heat of combustion (H,) - 

The ratio of laser intensity after passing through the 
flame and the initial laser intensity is known as the 
transmittance through the flame: 


ag “Ayse 
ae! lola eae Cea 
Ag be 


Under the present approximations, the transmittance is related 
to the volume fraction and the refractive index of scot along 
the radiation path. The refractive index of soot is assumed 
to be independent of temperature and fuel type. There is 
controversy in the literature concerning the temperature 
dependence of refractive index of soot ([Sivathanu, et al. 
(1991)]. However, variation with temperature and fuel type is 
not seen for a majority of the measurements. MTherefore, a 
constant value (1.55- 0.561) of refractive index given by 
Dalzell and Sarofim (1969) is used in the present study. With 
this value; (Ky g1sessetwrat 4r39. The laser absorption 
measurements are used to calculate 1,, using Eqn. 2-2. 
Radiation intensities leaving the flame at the two 
wavelengths (A = 0.9 and 1.0 pm) are related to the 
distribution of soot volume fractions and temperatures along 
the narrow radiation path selected by the long receiving tube 


by: 
Ty (S) 

,-f Ty dty= (lots). ty 12) (2-4) 

wh 

ere the local transmittance 1,(s) is defined by Eqn. 2-2 and 
I,, is the black body intensity at the local temperature: 
2nce 
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the two emission wavelengths defined by Eqn. 2-3 depends on 
the soot volume fraction distribution along the path weighted 
by the local Planck's function I,,- The K, is 5.33 at 0.9 um 
and 5.59 at 1.0 um based on the refractive index for soot. 

Equivalent radiation temperature (T.) and emission 
transmittance (1,,) for the radiation path are defined as shown 
in Eqn. 2-4. Using the measurements of intensity at the two 
wavelengths, Eqns. 2-4 are solved for T, and 17,,, assuming that 
the reciprocal wavelength dependence oOf.aEGn. 42-3 is 
applicable to 1T,,. For an isothermal path, T, is equal to the 
actual temperature and 7T,, is identical to 71,, after correcting 
for the wavelength dependence of the absorption coefficient. 
For a non-isothermal path, Sivathanu et al. (1991) found that 
the-T. is close to the peak soot temperature along the path. 
Ifa farge portion of the soot particles along the path are at 
temperatures’ close’ tol. ,Wthen 7,) is similar to 7,,. On/the 
other hand, if 1,. is much larger than T,,, a large portion of 
the soot particles are at relatively low temperatures. 

2.2.4 Operating Conditions 

The operating conditions for the present tests are 
summarized in Table 2-2. The general flame shapes were 
typical of pool fires with the flame width equal to that of 
burner at the lowest regions of the flame and then necking in 
toward the pool center at some distance above the burner. The 
fuels of similar composition, such as fuels classified as 
paraffins, aromatics or alcohols, generally have comparable 
flame shapes. The alcohol fuels generally necked in just 
above the burner lip and the parafin and aromatic fuels had a 
larger distance before necking begins. All flames exhibited 
large scale pulsations characteristic of buoyancy dominated 


flows (Pagni, 1989). The parafin and aromatic flames formed 
soot particles at all scales, which is characterized by yellow 
luminous radiation. The alcohols contained no measurable 


soot, though the ethanol and hexanol flames did exhibit faint 
visible luminosity. Methanol flames exhibit only radiation 
from the combustion gas bands and burned with a bright blue 
color. Toluene (an aromatic), heptane (a parafin), and 
methanol (an alcohol) were selected as representative fuels 
for more detailed study. These three fuels represent a wide 
range in sooting tendency. 


2.3 RESULTS AND DISCUSSION 

2.3.1 Radiative Properties 

Measurements to determine the radiative heat loss 
fraction for a variety of fuels yielded extensive information 
on the distribution of radiative heat fluxes to the 
surroundings. Measurements of radiant heat flux incident on 
the plane of the burner lip in the radial direction and 
incident in the axial direction at a fixed distance from the 
burner center are presented (see Fig. 2-4). As expected, the 
overall magnitude of the radiant heat fluxes varied with both 
fuel and scale. However, the trend of the data was remarkably 
similar for all fuels and scale. The general trends are shown 
for both luminous and non-luminous flames in the following 
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Table 2-2 Test Conditions for Liquid Fueled Pool Flames 


figures. Fig. 2-8 shows the distribution of axial heat fluxes 
for heptane at three different scales. Fig. 2-9 shows the 
axial fluxes of methanol for at the same three scales. The 
non-dimensional axial radiative heat fluxes, (UAB EON, are 
shown as a function of axial distance normalized by the flame 
height (x/H,;). The maximum flux was found at a relatively low 
axial position at all scales. The data collapses until x/H, 
= 0.8 beyond which scatter is seen. This trend that shows 
that radiative heat flux did correlate well with total heat 
release of the flame for luminous and non-luminous flames 
until the magnitudes of the heat flux data are low. Hamins et 
al. (1991) found that measurements of total radiative heat 
loss based on single location measurements were most accurate 
when made in the range of x/H, = 0.2 to 0.6. This corresponds 
to the location where the heat flux values at different scales 
collapse when normalized using this method. 

The radial heat fluxes for heptane and methanol at 
several scales are shown in Fig. 2-10 and Fig. 2-11. The 
radial radiative heat flux is normalized by the square of the 
distance from the burner center divided by the total chemical 
energy of the flame (q"*(r-D/2)? ) / Qa, (r-D/2) is used as the 
length scale for the radiative flux "Since m(r-D/2)* represents 
the area on the radial plane to which the flame radiates. The 
radial distance from the pool center is normalized by the 
flame height (r/H,). The heptane data collapse along a single 
curve for all flames considered until the heat flux magnitudes 
become low. The methanol radiative heat flux do not collapse 
well. However, the magnitudes of the flux are low in this 
region. 

The radiative heat loss fractions, X, (the fraction of 
total chemical energy of flame lost as radiation) were found 
to be nearly constant for a particular fuel, and relatively 
independent of burner size (Table 2-2). This property has 
been seen in other highly-buoyant flames (deRis, 1979). 
Furthermore, each classification of fuel, i.e. parafin, 
aromatic or alcohol, also tend to have a constant radiative 
heat loss fraction, regardless of fuel or scale. The parafin 
and aromatic fuels generally lose approximately 30% of their 
energy as radiation while the alcohols lose approximately 20%. 
Flames burning an azeotropic mixture of toluene and ethanol 
had an X, between the two groups. 

The relatively constant radiative heat loss fraction for 
the various fuels was quite surprising in light of the 
tremendous differences in sooting tendency of the hydrocarbon 
fuels. For example, the soot production for toluene and 
heptane varied by an order of magnitude (Gore et al. 1991). 
Soot continuum radiation contributes much to the radiant 
output of flames, but the radiant energy released by the 
flames and the burning rates of the fuels vary by less than 
50% at any one scale. If all the soot present in the flame 
was contributing the flame emission, the total radiant output 
and the fuel burning rates should also vary by approximately 
an order of magnitude for optically thin flames. This will 
further investigated by path integrated measurements of 
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radiation temperature and transmittance. 

2.3.2 Effective Temperatures and Transmittances 

Measurements of radiation temperature and effective 
absorptive and emissive transmittance were made for heptane 
and toluene flames from the 7.1 cm burner. The mean emission 
temperatures for the two flames are plotted as a function of 
normalized distance from the fuel surface in Fig. 2-12. As 
noted in section 2.2.3, these temperatures are close to the 
peak temperatures across the entire diametric radiation path. 
For approximately half the flame height, the temperatures in 
the toluene flame are lower by approximately 200 K than those 
of the heptane flame. The adiabatic flame temperature 
(assuming complete combustion) for toluene is approximately 
150 K higher than that for heptane. The observed temperatures 
indicate incomplete and non-adiabatic combustion. Since the 
radiative heat loss fraction for the two fuels are close, the 
differences in the temperatures can not be attributed to the 
differences in radiative cooling but are consequence of the 
reduced combustion efficiency of toluene as compared to 
heptane. 

Eig 2-13 shows the absorption and emission 
transmittances for the two flames plotted as a function of 
normalized axial distance. The measurements of 1,, show that 
the toluene flames are optically thick with a transmittance of 
less than 5% at certain locations and generally less than 30% 
at all locations. Based on both 1,, and 1,,, the heptane 
flames are optically thin with transmittances greater than 70% 
for all locations. The differences between 1,, and T,, for the 
heptane flames are relatively small suggesting that large 
portions of the soot particles in this flame contribute to the 
radiation intensity. In contrast, the differences between 1,, 
and T,, for the toluene flame vary by a factor of 4 to 10. 
This observation suggests that a large portion of the soot 
particles in the toluene flames are at relatively low 
temperatures and do not contribute significantly to the 
radiation intensity. Sivathanu et al. (1990) found that a 
large percentage of soot present in strongly radiating 
acetylene/air jets was relatively "cold" and not contributing 
to flame emission. The existence of large quantities of cold 
soot in the toluene flames helps explain the relatively small 
increase in X, and heat feedback for fuels of different 
sooting tendencies. 

2.3.3 Flame Heights 

Visible flame heights for the test flames in the present 
study are plotted in Fig. 2-14 and tabulated in Table 2-2. 
Flame heights measurements from the 4.6 cm and 60 cm burners 
were not made in this’ study. The results are plotted 
according to the scaling relationships proposed by Zukoski et 
al. (1984) with H,/D as a function of normalized heat release 
rate (Q,.,/a/*). Data for a variety of buoyant diffusion 
flames from gaseous, liquid and solid fuels (Thomas, 1962; 
Steward, 1970; Terai and Nitta, 1975; McCaffrey, 1979; You and 
Faeth, 1979; Zukoski et al., 1984; and Sivathanu, 1990) are 
plotted as solid lines in Fig. 2-14. 
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The flame heights for methanol are characterized very 
well by this normalization and agree with the data from other 
investigators. The flame heights from heptane measured in 
this study are in the lower portions of the correlation for D= 
7.1 and 30 cm. The flame height for heptane at D = 100 cm and 
the normalized toluene flame heights at all scales do not 
correlate well as a function of normalized heat release when 
ideal heat of combustion is considered. This may be due to 
the low combustion efficiency of toluene (0.76 according to 
Mulholland et. al., 1988) which is not taken into account by 
the correlation. If the toluene flame heights are corrected 
for combustion efficiency (Q,. in Fig. 2-14), agreement with 
correlations from the literature improves. Additionally, 
blockage of the upper regions of the toluene flame by soot 
release (especially at D = 100 cm) made accurate flame height 
measurements difficult. 

2.3.4 Fuel Burning Rates 

Fuel burning fluxes for many of the present test flames 
are shown in Fig. 2-15 and tabulated in Table 2-2. The fuel 
burning rates have been plotted as a function of the burner 
diameter. As described previously, a constant fuel height was 
maintained in the burners during measurements. 

Several researchers have compiled fuel burning rate data 
versus diameter for a variety of fuels (e.g. Barbasukas, 1983; 
Hamins et al., 1992). Burning rate data of Corlett and Fu 
(1966) for methanol and Koseki and Yumoto (1988) for heptane 
are included in Fig. 2-15 for reference. Corlett and Fu 
(1966) maintained a constant fuel level relative to the burner 
rim. An increase of mass burning flux at the smaller 
diameters (D > 3 cm) has been found by other investigations 
(e.g. Blinov and Khudiakov, 1957). This may be due to the 
enhanced effect of conduction from the burner material to the 
fuel or laminar convection at small diameters. This effect 
was not studied in the present investigation. At larger 
diameters, the alcohol fuels show very little variation in 
burning flux rate with diameter. 

The parafin and aromatic fuels show an increased mass 
burning flux with diameter. This was attributed the increased 
importance of radiative heat feedback to the fuel surface as 
the pool diameter increases (Blinov and Khudiakov 1957, 
Hottel, 1958). Koseki and Yumoto (1988) allowed the fuel 
level relative to the burner rim to decrease with time. This 
may explain the systematic difference between their burning 
rate data and the present study. 

2.3.5 Overall Energy Distribution 

Table 2-3 shows the energy distribution of the 30 cm pool 


fires of methanol, heptane and toluene. The ideal heat 
release (Q,.,) and the total feedback (Q,) are based 3 the 
burning rate of the fuel. The energy ‘Poss term (Q, is 


estimated from the increased sensible energy of the euesy and 
the increased temperature of the cooling water. The sensible 
energy of the fuel was estimated from the liquid fuel 
temperature distribution. 
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2.3.6 Spectral Fuel Vapor Transmittance and Flame 
Emission 

The spectral content of flame emission at x/D = 1.0 for 
heptane and toluene are presented in Fig. 2-16. The greater 
amount of soot in toluene accounts for the generally larger 
emission signal as compared to heptane. The spike at 4500 nm 
is due to CO, and is quite evident in both flames. In 
heptane, the emission at 4500 nm is quite large when compared 
to soot emission which is evident across the spectrum. This 
indicates that gas band radiation in heptane should be quite 
significant. In contrast, the CO, emission at 4500 nm for 
toluene, though apparent, is not overwhelming compared to the 
soot emission across the spectrum. 

Fig. 2-17 presents the fuel vapor transmittance for 
heptane and toluene. Significant drops in transmittance are 
seen at wavelengths of approximately 3400 nm for both heptane 
and toluene and 4500 nm =*{for heptane. The drop in 
transmittance at 4500 nm for toluene may have been missed due 
to the large (200 nm) measurement interval as opposed to the 
heptane measurements. The 4500 nm line corresponds to the 
presence of CO, in the flame. The drop in fuel vapor 
transmittance does not correspond with regions of high flame 
emission (Fig. 2-16) indicating that fuel vapor above the fuel 
surface is not a important absorber of incident radiation. 

2.3.7 Fuel Surface Reflection 

A possible component of the radiative feedback process in 
pool fires could be the reflection of the incident radiation 
off the fuel surface. Fig. 2-18 shows measurements of thermal 
radiation off the surface of water, toluene, methanol, and 
heptane at various incident angles. At very low incident 
angles, a large percentage of the incident energy is reflected 
off the fuel surface and is not involved in the vaporization 
of the fuel. 


2.4 Summary and Conclusions 

13 The total radiative heat output for a variety of 
liquid fuels was obtained using a multi-point measurement 
technique. Several alcohol fuels were considered at different 
scales and a radiative heat loss fraction of approximately 0.2 
was found. A variety of parafin and aromatic fuels, with 
vastly different sooting tendencies, were also tested at 
different scales and a radiative heat loss fraction of 
approximately 0.3 was found over the range investigated. 


Ze The radiative heat flux distributions for parafin 
and aromatic fuels correlated fairly well as a function of 
heat release rates and flame length in the axial and radial 
directions. Radiative heat flux for methanol correlated well 
in the axial direction, but showed a large amount of scatter 
in the radial direction possibly due to the small levels of 
Signal. 


3% The effective radiation temperatures and 
transmittances for heavily sooting toluene and moderately 
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FLAME EMISSION (W/m?—sr—pm) 
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Fig. 2-16 Flame emission at various wavelengths for toluene and 
heptane flames 
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Fig. 2-17 Fuel vapor transmittance at various wavelengths for 
toluene and heptane 
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sooting heptane indicate that toluene contains a large volume 
of "cold" soot particles which are not contributing to 
emission. This may explain the similarity of X, and heat 
feedback for fuels of differing sooting tendency. 


4. The visible flame heights measured during the 
present investigation satisfied the correlations of Zukoski et 
al. (1984) for buoyant diffusion flames. Flame heights for 
methanol correlated very well at all scales. At the larger 
scales, flame heights for heptane and toluene did not 
correlate well as a function of idealized heat release, 
possibly due to low combustion efficiency. Accurate 
measurement of toluene flame heights was difficult at larger 
scales due to blockage by cold soot particles. 


55 Measurements of flame emission and fuel vapor 
transmittance in toluene and heptane indicate that fuel vapor 
absorption of radiation is not significant in these flames. 
Flame emission measurements show that gas band radiation in 
heptane flames should prove to be important, while in toluene 
flames, continuum soot radiation is dominant. Reflection of 
incident infrared radiation (A = 1.0 - 3.0 um) is substantial 
at low incident angles. 
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CHAPTER III 
RADIATIVE FEEDBACK TO THE FUEL SURFACE 


3.1 Introduction 

A major goal of this investigation is to improve the 
understanding of the role of radiation in providing energy for 
fuel vaporization in liquid pool flames. Measurements of 
radiation heat flux incident on the fuel surface in different 
directions at various radial positions are necessary for this 
purpose. In the following, existing measurements of 
hemispherical radiative heat feedback and the techniques used 
for these are reviewed. Measurements of directional heat 
feedback to the fuel surface have not been reported in the 
past. Therefore, a new technique was developed for obtaining 
these. The technique and its application to the pool fires is 
described in this chapter and the resulting measurements are 
reported and discussed. Fuels of various sooting tendencies 
were used and differences in radiative feedback 
characteristics are examined. 

Blinov and Khudiakov (1957) monitored burning rates for 
a variety of industrial fuels in pools consisting of 
concentric rings. Total heat transfer to the fuel within each 
ring can be found based on the local fuel burning rates within 
each ring. Blinov and Khudiakov (1957) found the highest heat 
transfer at the pool center, decreasing somewhat along the 
radius, and increasing at the pool edge. Akita and Yumoto 
(1965) used a multi-ring pool to obtain local measurements of 
burning rates in methanol. MThis study found qualitatively 
different results than Blinov and Khudiakov (1961) for their 
methanol pools. Hamins (1992) found results similar to Blinov 
and Khudiakov (1957) for several fuels, including methanol. 
Modak and Croce (1977) measured the local fuel surface 
regression rates for various sizes of square PMMA fuel beds 
(.025 - 1.22 m (Side)). Radiative feedback was found from a 
surface energy balance. Convective heat feedback was assumed 
to be the primary heat transfer mechanism for the smaller 
diameter pools and taken as a constant. Heat feedback was 
shown to be dominated by thermal radiation at larger scales. 


A few investigators have reported direct measurements of 
heat feedback at the fuel surface in order to determine the 
importance of radiation and convection components. Corlett 
and Fu (1965) measured the fuel vaporization rate in an 
insulated teflon cup containing fuel placed in 30 cm pools. 
Their calorimetric measurements indicated that radiation was 
an important component in non-sooting methanol fires. Alger 
et al. (1979) used Gardon and transpiration radiation gauges 
to measure total and radiative heat feedback in large (3 m) 
methanol pool fires. Shinotake et al. (1985) used gardon 
gauges with differing surface emissivities to separate 
convective and radiative components. Measurements at four 
locations in heptane pool fires of 100, 70, and 30 cm were 
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made. Results indicated that radiation contributed 80% of the 
total feedback for the 100 cm and 80 cm pools and 70% of the 
total feedback for the 30cm pool. Yumoto (1971) also used two 
color Gardon gauges in intermediate sized hexane pools. This 
study showed slightly larger contributions of radiation to the 
total heat feedback. Changes in surface emissivity of the 
gauges by fuel condensation and soot deposition are some the 
problems associated with this technique. 


3.2 Experimental Methods 

3.2.1 Experimental Apparatus 

The details of the pool burner and fuel supply system are 
described in Section 2.2. Radiative feedback measurements 
were performed in the 30, 60 and 100 cm diameter burners. A 
constant fuel level (5 mm below burner rim) was maintained. 
The apparatus was allowed an adequate warm-up period before 
data were acquired. The warm time was determined by a trial 
and error procedure to obtain steady burning rates and flame 
heights. The feedback radiometer was mounted on a traverse 
system that allowed for linear motion in two directions and 
angular motion about the radiometer axis. The traverse system 
enabled measurements to be taken at locations over most of the 
fuel surface. Measurements of radiative feedback in the 60 cm 
burner for heptane and toluene were made in the center of the 
burner due to traversing difficulties caused by the large heat 
release of these flames. Feedback measurements in the 100 cm 
burner were made for methanol only. 

3.2.2 Instrumentation 

The radiometer was specifically designed for use in the 
hostile environment near the fuel surface. The radiative 
component of heat feedback was measured independently avoiding 
the need for approximate models for deducing this quantity. 
Convective heat transfer was eliminated by placing the heat 
sensing element inside a gauge housing which is water cooled 
and has a purge flow of nitrogen over the sensing element and 
through its aperture. 

Fig. 3-1 shows a schematic of the gauge. The sensing 
element is recessed 11 mm below the top surface of the gauge. 
This restricts the gauge field of view to 14° (half peak). 
This design allows measurements of directional heat fluxes 
resolved to a relatively narrow angle. The methodology for 
determining total hemispherical radiative heat flux from 
directional measurements is fully described in section 3.3.3. 
The spectral response of the gauge was flat between 0.5 and 10 
um. A nitrogen purge (30 cc/min) was used to prevent soot and 
fuel vapor deposition on the sensing element. Care was taken 
to determine that this purge had minimal effect on the flame. 
Additionally, considerable cooling the gauge using tap or warm 
water was required to maintain the gauge at a relatively 
constant temperature. Increased gauge temperature changed the 
baseline zero. Fig. 3-2 shows the arrangement used for 
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traversing the gauge, changing its view direction and 
delivering the purge and cooling water flows. Cooling water 
and the nitrogen purge were delivered to the gauge through 
tubing which ran through the fuel residing in the burner. This 
design allows for minimum intrusion of the gauge on the flame. 

The signal from the radiometer was sent to a digital 
amplifier (gain = 500) and then sampled by a data acquisition 
computer. The sampling rate was 30 Hz. The signal was 
sampled for 90 seconds at each location. Effects on the 
Signal by gauge heating were found by measurement of the 
baseline signal. This was accomplished by placing a water 
cooled cap over the radiometer (with nitrogen purge 
maintained) with the flame burning but the flame radiation 
blocked from the gauge. An uncertainty of 10% for the 
intensity measurements is estimated from repeated trials. 

Radiometer Calibration 

Calibration consisted of comparison of the voltage 
generated by the radiative feedback radiometer with that 
generated by a standard wide angle radiometer (150° field of 
view) in response to a known blackbody flux. Measurements of 
incident intensity from the black body were made as the 
feedback radiometer was rotated about the center of the top 
surface. These values were integrated using the technique 
described in section 3.3 and compared with the standard 
radiometer to obtain the calibration factor. Nitrogen purge 
and the appropriate water cooling were used during 
calibration. 


3.2.3 Data Analysis 

Due to the narrow field of view, the feedback radiometer 
measured incident directional intensity at a given location on 
the fuel surface. The hemispherical radiant heat flux to a 
point is found by integrating all directional intensities 
intersecting a hemispherical surface placed over that point. 
Fig. 3-3 (following Siegal and Howell (1981)) indicates the 
geometry used in this integration. Intensity measurements are 
made at several positions about the surface of the hemisphere. 
These values are then integrated over the hemisphere using the 
trapezoidal rule utilizing (Siegal and Howell, 1981): 


Oraa~/ [r(8,)] sin(8) cos (8) dbd (3-1) 


The quantity sin © de d¢@ in Eqn. (3-1) describes the surface 
area on the unit hemisphere that the intensity in direction 
(8,¢) intersects (Fig. 3-3). 

Axial symmetry allowed the integration to be separated 
into two parts. The limits of integration over two quadrants 
of the hemisphere are defined by: 


m/2 and 7/2 < @ < 31/2 
Thee candi] 2e<yhr <r / 2 (3-2) 
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where © = 0, ¢= O defines a ray parallel to the fuel surface 
towards the pool center and ® = 90, ¢@ = O defines a ray normal 
to the fuel surface. 

3.3.4 Operating Conditions 

Radiative feedback measurements were made for three 
fuels: toluene, heptane, and methanol. The test conditions 
are summarized in Table 3.1. Intensity measurements were made 
at a minimum of 8 radial positions across the pool surface. 
At most radial locations, angular distributions of a minimum 
of 11 pairs of angles (@,¢) were found (before assuming axial 
symmetry). Where direct measurements were not possible, 
linear interpolation was used to fit an intensity value for 
that position. 

Tap water was used for cooling the feedback gauge during 
the heptane and methanol tests. In toluene, vapor 
condensation on the gauge surface occurred when cooling with 
tap water was utilized. Therefore, warm water (approximately 
60°C) pumped from a water bath was used for cooling the gauge 
during these tests. As the water lines were routed through 
the liquid fuel, the temperature of the cooling water affected 
the fuel burning rate. This effect is noted in Table 3-1. 
Signal baseline changes were measured during each data 
collection session. The baseline did change with gauge angle 
but proved to be fairly insensitive to radial position along 
the fuel surface. 


3.3 Results and Discussion 

3.3.1 Radiative Feedback in the 30 cm Burner 

Figure 3-4 shows directional radiative intensity as a 
function of distance for the methanol pool fire for 
representative directions. The angular variation is shown for 


angles looking away from pool center (¢@ = 0). Figure 3-5 
shows a similar figure for the heptane pool fires with the 
same representative directions as shown in Fig. 3-4. The 


burner center is located at r = 0 cm. It was found for all 
fuels and radial positions, the peak intensity occurred at 6 
= nm/2 when the radiometer was looking away from the pool 
center (@ = 0). When @ = 7m, the intensities initially 
increased with distance from the pool center. This is due the 
increased radiation path length through the flame. Figure 3-6 
shows the intensity distribution for the heptane flame at r = 
0 cm (pool center) for various angles of © at ¢= O and 7. 
Figure 3-7 shows a similar plot for heptane at r = 8 cm. 
Significant intensities were measured from all directions. 

The intensities from the heptane flame are twice as large 
as those from the methanol flame for a given angle and radial 
position. This is due to contribution of continuum radiation 
from soot in the heptane flame, which is not present in the 
methanol flame. Toluene showed similar magnitudes of 
intensity as the heptane flames for directions normal to the 
fuel surface but received greater contributions from angles 
off the normal to the fuel surface. 
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"estimated based on complete combustion to stable gaseous 
products 


Table 3-1. Operating Conditions for Radiative Feedback 
Measurements 
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The incident intensity values are integrated over the 
hemispherical surface using the technique described in section 


3.3. Figure 3-8 shows the incident radiative heat flux for 
toluene, heptane and methanol as a function of radial 
position. The radiative heat flux profiles for toluene and 


heptane are fairly flat profiles across the fuel surface. In 
heptane, the radiant heat flux drops slightly only at the edge 
of the pool. However, methanol flames have a large reduction 
in the radiant heat flux from the pool center (r = 0 cm) to 
the pool edge. This result is evident in the flame shapes of 
the various flames. Toluene and heptane flames are quite 
broad at the base and neck towards the pool center 
approximately 5 cm above the burner. This shape provides a 
consistent radiation source across the majority of the fuel 
surface. The methanol flame swept across the fuel surface 
just above the burner lip , leaving a thin flame over the fuel 
at large radii. The flame begins to neck upwards only near 
the center of the burner. 

Hamins (1992) measured burning rates in 30 cm concentric 
ring pool flames for methanol, heptane and toluene. The 
steady-state heat flux incident at the fuel surface in each 
ring was: 


qm! (H,po+C, Gee) “ye; (3-3) 


where i is the ring number, q," is the average energy flux 
(kW/m?) incident on the i‘ annular ring, m", is the average 
mass burning flux (g/m*-s) in the it? annular’ ring, Hy, is the 
heat of vaporization, T,; is the pool surface temperature 
assumed to be equal to the boiling temperature, T, is the 
ambient temperature, and q." is the sum of the energy lost to 
the surroundings and the increased sensible energy of the fuel 
resident in burner. The sensible energy increase was 
calculated from the temperature gradients measured in the 
liquid fuel. The energy lost to the surroundings was found by 
measuring the temperature change in the cooling water 
reservoir. The loss term q," was measured to be 1%, 5% and 8% 
of the incident energy (Q) for the methanol, heptane, and 
toluene fires respectively. 

Figure 3-9 compares the results from the radial variation 
of radiation feedback measurements normalized with the burning 
rate measurements of Hamins (1992). The radiative feedback 
results have been corrected for fuel burning rate differences 
caused by the presence of the radiometer. The experimental 
uncertainty of the normalization presented in Fig. 3-9 is 
estimated to be 20%, mainly due to the effects of different 
burning rates in the tests. Figure 3-9 shows that the heat 
feedback is dominated by radiation across the entire burner 
for toluene and heptane fires. Methanol also has a 
significant radiation component, even without the presence of 
soot. At the pool center, heat feedback in methanol fires is 
also dominated by radiation. Convection played a significant 


49 


quanjo) puke suejday ‘joueyoW 10} UOIISOd [eIpel JO UONOUNJ & Se YORqpaaj yeoY sANeIpPeEY g-¢ ‘SIq 


aDdaq 100d Wd “UOT}ISO” [eIpey YALNAO 1OOd 
Gl cl 6 9 € 0 


ouonjoy ov 
ouridsayy so 
jourma,—@ 


W/ MAY “XN]J OH OANLIPeY 


50 


(%0c + Ayureyigoun) 


guaN{O} pue suR}doy ‘JouRYyoUI 10J UOTISOd [eIped JO UOIJOUNJ B Se YORGpaaj Je9Y VANVLIPeI PIZI[BUION 6-¢ “Sly 


dda TOOd 
Gl 


col 


Wd ‘UOTJISOg [eIpeYy 
6 9 


YAaLNAD TOOd 
€ 0 


00 


ouanjoy, =v 

ounjdoyHH sod 

joueya, se 

C0 

) 
& 
~~ 
=) 
2 


Ol 


Cl 


5] 


Incident Radiative Flux (kw/m* ) 


Fig. 3-10 Radiative heat feedback as a function of radial 
position for methanol in 60 and 100 cm burners 
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role near the pool edge, where radiation contributed 35% of 
the total heat feedback. 

3.3.2 Radiative Feedback in 60 and 100 cm Burners 

Measurements of incident radiative feedback to the fuel 
surface for methanol (60 and 100 cm) are presented in Fig. 3- 
10 as a function of radial position. As was found for 
methanol in the 30 cm burner, the incident radiation decreases 
as the distance from the pool center increases. 

Measurements of incident radiation on the fuel surface at 
r/D = 0.0 for toluene found 47.6 kW/m’. This value is 
considerably higher than the average amount of energy 
necessary for vaporization of the fuel (28 kW/m* based on 
burning rate). 


3.4 Conclusions 

Radiative heat feedback was measured using a narrow- 
angle, nitrogen purged radiometer which was placed just above 
the fuel surface. Considerable care was taken to prevent fuel 
vapor condensation and soot deposition on the sensing element 
of the radiometer. The angular distribution of radiative 
feedback was found for positions across the radius of a 30 cm 
pool for fuels of differing flame luminosity and shape. The 
main conclusions from the radiative feedback measurements are 
as follows: 

1. Considerable contributions to radiant feedback flux 
are found from all directions. Multi-ray calculations of 
radiative heat feedback must consider the intensity 
contributions from different directions. 


2h Radiation was found to be a significant, heat 
transfer mechanism in medium and large scale pool flames of 
vastly different luminosity. The radiative feedback 


contributed almost the entire heat transfer to the fuel 
surface in toluene and heptane fires. Radiative heat feedback 
was also found to be large near the center of methanol pool 
fires, but the radiative contribution diminished considerably 


near the pool edge. 


The measurements discussed in this chapter show the 
importance of radiative heat transfer for fuel vaporization in 
medium-sized pool flames. Measurements are available for the 
evaluation of predictions of radiative feedback using local 
flame emission and absorption characteristics. Such an 
evaluation is undertaken in Chapter 5. 
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| CHAPTER IV 
LOCAL EMISSION/ABSORPTION PROPERTIES 


4.1 Introduction : 

Emission and absorption properties of liquid pool flames 
are of interest in fire safety applications due to their 
strong influence on the radiant heat loading and radiant heat 
feedback. Knowledge of heat feedback from flames to fuel 
surface is also relevant in other applications such as spray 
flames and solid rocket motors. The objective of the present 
work was to obtain and interpret simultaneous’ local 
measurements of spectral emission and absorption of radiation 
in liquid fueled pool flames. Toluene was selected as the fuel 
based on its frequent use as an industrial solvent and because 
of its strong sooting tendency which is common to most 
practical hazardous fires. 

Experimental studies of pool flames date back to the work 
of Blinov and Khudiakov (1957) and Hottel (1958). These 
authors recognized the turbulent nature of pool fires at 
relatively small pool sizes. The importance of radiative heat 
transfer in determining the burning rates at relatively large 
scales was recognized (Hottel, (1958) and Burgess and 
Hertzberg, (1974)). The early models utilized an average flame 
emissivity, a shape factor and a constant flame temperature in 
the calculation of radiation heat flux (see Burgess and 
Hertzberg, (1974) for example). All of these parameters were 
obtained from experimental data regarding burning rates in 
conjunction with a homogeneous, gray flame assumption. Modak 
(1977) developed a method to obtain the individual shape 
factors for various radiation paths through the flames and 
treated the absorption coefficient and the effective radiation 
temperatures as properties of a complete radiation path. The 
effective temperature was obtained from Schmidt measurements 
and the local gray absorption coefficients were obtained from 
blackbody absorption data. Modak (1981) upgraded this method 
to include local measurements of mean temperatures and mean 
absorption coefficients obtained from path integrated 
blackbody emission and absorption data and subsequent 
deconvolution. Modak's models (1977, 1981) and other more 
recent global models (Shorki and Beyler, (1989) and Hamins et 
al. (1992)) rely on the assumption that the effects of 
fluctuations in temperature and absorption coefficients on the 
radiation heat flux (termed turbulence-radiation interactions- 
see Faeth et al. (1989) for example) can be treated using 
effective mean radiation properties. Orloff (1981) used 
Markstein's (1981) scanning radiometer measurements for 
obtaining effective temperature and absorption coefficients as 
a function of distance from the fuel surface. Variations in 
radiation intensity from different radial positions were 
assumed to originate from different view angles and flame 
shapes. The mean flame shapes were obtained by averaging 
instantaneous photographs similar to the work of Modak (1977, 
1981). 
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Markstein (1981) obtained instantaneous measurements of 
emission intensity and absorption coefficients for arbitrary 
paths in pool fires. These data were processed to obtain 
moments of intermittency and spatial correlations of radiation 
intensity fluctuations. Markstein (1981) suggested that such 
information regarding turbulent properties was necessary for 
improved models of radiation from pool flames. Effects of 
turbulent fluctuations on radiative heat losses from gaseous 
jet and pool diffusion flames have been well-documented (see 
Kounalakis et al. (1988a,b), Sivathanu et al. (1990), Gore et 
al. (1991) and their references). 

Fischer et al. (1987) measured the fluctuating 
temperature and mean absorption coefficient in a 0.5 m ethanol 
pool fire. Radiation predictions utilizing the measured mean 
properties were consistently lower than the experimental data. 
Enhanced predictions of intensity were obtained § by 
incorporation of measured temperature fluctuations with the 
nonadiabatic fast chemistry approach. This further 
highlighted the need for including the effects of scalar 
fluctuations in radiation models. 

The present study is focussed on determining the scalar 
property distributions including turbulent fluctuations in 
toluene-fired pool flames. The effects on heat feedback are of 
particular interest since the radiative flux in this case may 
originate in the visibly persistent flame zone near the liquid 
surface. Recently, Bouhafid et al. (1989) have detected 
intermittency in the visibly persistent zone. However, these 
authors suggest that the characteristics of the fluctuations 
near the liquid surface are those of an oscillatory laminar 
flame rather than a fully turbulent flame. It is therefore of 
interest to measure mean and turbulent scalar property 
distributions over an entire liquid-fired pool flame and then 
to evaluate the effects of fluctuations on heat flux incident 
at various locations including the fuel surface. 
4.2 Experimental Methods 

4.2.1 Experimental Apparatus 

The details of the pool burners and fuel supply system 
are included in section 2.2. The emission/absorption 
measurements were made in the 7.1 cm burner. A constant level 
of fuel 3.5 mm below the burner rim was maintained in the pool 
burner. 

4.2.2 Instrumentation 

Emission and absorption measurements resolved to 6 mm 
were obtained throughout the flame using the three-line 
technique described by Gore and coworkers (Sivathanu and Gore, 
(1992), Sivathanu et al. (1991) and Sivathanu and Gore 
(1991)). Wideband emission intensities were measured at 0.9 
um + 50 nm halfwidth and 1.0 wm + 50 nm halfwidth and 
transmittance of a laser beam was measured at 0.632 um. 

The three-line probe was a slightly modified version of 
those used by Gore and coworkers (Sivathanu and Gore, (1992), 
Sivathanu et al. (1991) and Sivathanu and Gore (1991)) and is 
shown in Fig. 4-1. The local emission signals were collected 
by a 6 mm diameter tube which was slightly pressurized by a 1 
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cc/min flow of nitrogen. The light passed through the window 
at the exit and through a pair of beam splitters. The first 
beam splitter only served to send the 0.632 um laser beam into 
the optical probe in an opposite direction. The second beam 
splitter separated the emission beam into two parts which were 
then incident on the two photomultiplier tubes equipped with 
the wide band monochromatic filters. The direction of the 
laser beam was reversed in order to allow the operation of the 
PMTs at high sensitivity in order to extend the range of 
temperature estimates to approximately 850 K. The wide band 
filters (particularly the one at 0.9 um) have non-zero 
transmittance at 0.632 um leading to interference from the 
laser light when the PMTs were powered for high sensitivity. 
Reversing the direction of the laser beam avoided this 
interference. Since the laser transmittance data were not used 
in the interpretation of the emission signals, the reversed 
direction of the two light signals did not affect the results. 
The laser beam was collected by a second purged probe after 
passing through the probe volume. 

Where ever possible, non-intrusive optical probes are 
preferred in comparison to the present intrusive probes. 
However, existing non-intrusive methods are not effective in 
the present optically thick unsteady flames. Therefore, the 
effects of the intrusive probes and purge flow on the local 
properties and effects of the probe heat up on the signal were 
evaluated by Sivathanu et al. (1991). 

The local temperature and mean species concentrations 
were affected by less than 5% by the introduction of the probe 
as summarized in Sivathanu et al. (1991). The nitrogen purge 
flow was photographically observed to produce 1 mm thick dark 
zones on the surface of the probes. Measurements at 
representative locations with double the purge flow rate (2 
cc/min) showed less than 5% differences in the emission 
intensities and absorption transmittances providing an 
estimate of the errors. 

A stainless steel plate was used to cover the front end 
of the receiving probe to observe the effect of probe heat up 
with the flame in place on the signal at the end of each 
measurement. If the signal due to probe heat up exceeded 5 % 
of the local mean emission signal, the measurements were 
repeated with shorter data collection time but more 
repetitions. The maximum data collection times were less than 
5 minutes due to the above probe heat up criterion. 

The signals collected by the PMTs and the laser power 
meter were amplified and processed through eighth order 
antialiasing analog Butterworth filters set at 125 Hz before 
sampling by an A/D converter and a laboratory computer at 250 
Hz for post-processing. A minimum of 5000 data points 
(occasionally from different realizations due to probe heat up 
limitations) were used to obtain averages. The sampling 
frequency was determined by power spectral analyses of the 
signals from representative locations. The sensitivity of the 
optics and the detectors was calibrated using a commercial 
blackbody operated at 1273 K before and after each test. The 
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two probes serve as relatively cold background for each other. 


4.2.3 Data Analysis 

4.2.3.1 Emission/Absorption Data 

The laser transmittance signals were used to estimate the 
soot volume fractions using the Rayleigh limit for small 
particles. Past multiline measurements (Gore and Faeth, 
(1988)) have suggested that this approximation is adequate for 
the present flames. Following Tien and Lee (1982), the soot 
volume fractions were calculated as: 


f,,7-1n(I,/1,°) A/K,S (421) 


The argument of the natural log in Eqn. 4-1 is the measured 
transmittance. "S" is the separation distance between the 
probes. The separation distance was 10 mm. The constant K, in 
Eqn. 4-1 involves the complex refractive index of the soot 
particles (n-ik) at the measurement wavelength (Dalzell and 
Sarofim (1969)): 
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The complex refractive index of soot particles has been 
measured by several investigators and considerable discussion 
of the effects of temperature, fuel type and wavelength exists 
in the current literature (see Sivathanu et al. (1990) and 
Sivathanu et al. (1991) for example). Consistent with earlier 
work (see Sivathanu et al. (1990) and Sivathanu et al. 
(1991)), the values given by Dalzell and Sarofim (1969) were 
used resulting in K, = 4.89 at a wavelength of 0.632 um. The 
soot volume fractions obtained from Eqn. 4-1 are designated as 
f,, to indicate that these were obtained from an absorption 
measurement. 

The interference filters associated with the emission 
detectors have a relatively large spectral width (+50 nm at 
half peak) requiring the use of integration over the filter 
transmittance curve. The intensities measured by the emission 
detectors are related to the local properties as: 


K, 


A2 
[) a-t) Aad 
Si, eae a ri (4-3) 
2 

F,da 
a1 
where T, is the monochromatic transmittance of the material in 
the probe volume over the probe length L, F, is the 
characteristic transmittance curve of the interference filter 
with lower and upper cutoff wavelengths of 41 and A2. I,, is 
Planck's function at the probe volume temperature T: 


2hc? 


7 Pn BET = 
>» 13 (exp (hc/AK,T) -1) Sela 


For the emission measurements at 0.9 and 1.0 um, the soot 
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particles in the probe volumes were the dominant contributors 
to the emission intensity and hence 1,. Therefore for these 
two wavelengths the transmittance in Eqn. 4-2 was equated to 
the transmittance due to soot particles: 


Ky fyeS 
ey Ors 
Ty g7e 


where K, = 5.58 for A = 0.9 wm and K, = 5.47 for A = 1.0 um 
are constants dependent on the refractive index of soot [20]. 
Since h, c and k, are fundamental constants and K, for A=0.9 
um and A=1.0 um are known, Eqns. 4-2- 4-4 for these two 
wavelengths were solved for the two unknowns - temperature T 
and soot volume fractions based on emission - f,,. Due to the 
integral nature of Eqn. 4-2, an iterative solution must be 
obtained. 

If the soot particles in the probe volume had a uniform 
temperature then the measurements soot volume fractions based 
on absorption and those based on emission would be identical. 
If the temperatures of the soot particles were different, then 
f,. represents a Planck function weighted average of all the 


particles. Therefore, f,. is biased towards the hottest soot 


particles in the probe ‘volume WHLLE Sr. 1s a measure, of all 
soot particles. 

As summarized in Sivathanu et al. (1991), the absolute 
uncertainties of the present data depend on the values of 
refractive indices for the soot particles. For a fixed set of 
refractive indices, the uncertainties in the mean temperatures 
are less than 40 K for f,, above 0.1 ppm and increase to 80 K 
for f,, = 0.01 ppm. The uncertainties inet edatagare tessethan 
10% at 10 ppm and increase to 40% at 0.01 ppm due to 
digitization errors. The uncertainties in the emission soot 
volume fractions are dominated by the choice of the refractive 
indices and are estimated to be 30 % at 1400 K and f,, = 1 ppm. 
The uncertainties in the measurements of RMS values were 
estimated by repeating the measurements at several 
representative locations. The uncertainties in RMS 
temperatures varied between 70 and 150 K, the uncertainties in 
RMS measurements of f,, varied between 25% at 10 ppm to 100 % 
at 0.01 ppm. The uncertainties in the RMS T eOocemaLe 
estimated to be less than 40% based on repeated measurements 
at representative locations. 

At several axial locations in the flame, monochromatic 
radiation intensity measurements were completed for 
progressively larger probe volume lengths. These measurements 
were accomplished by increasing the distance between the 
probes, at intervals of 10 mm, from the initial distance of 10 
mm until the entire flame width was within the probes. Mean 
and RMS monochromatic radiation intensities are then obtained 
as function of radiation path length. This procedure allows 
an examination of the individual terms in the equation of 
radiative transfer with emphasis on turbulence-radiation 
interactions. 

4.2.3.2 Turbulence-Radiation Interactions 


(4-5) 
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The radiation intensity and heat flux to surrounding 
locations leaving the flame depend on the contributions from 
all the local segments. The local measurements were obtained 
with a resolution of 10 mm. Whether this spatial resolution 
is adequate for the estimate of radiation intensity depends on 
the local turbulence length scales. 

If average temperature and soot volume fractions are used 
in the calculation of the emission intensity using Eqn. 2, the 
effects of nonlinear dependence of Planck's function on 
temperature and the nonlinear dependence of transmittance on 
soot volume fractions can cause significant errors. 
Furthermore, cross correlations between 1, and I,, are 
neglected by this approach. The relative importance of these 
nonlinearities needs to be examined for pool fires. 

A third consideration of the turbulent properties 
involves the spatial correlations between the scalar 
properties at adjoining and separated segments. Since 
instantaneous, multi-point data are not available, only 
inferences concerning the adequacy of the present single point 
data can be drawn. 

4.2.3.3 Frequency Analysis 

In order to examine these mechanisms, the instantaneous 
equation of radiative transfer is averaged to obtain equations 
for mean and mean square intensities at various points. The 
results of these equations were compared with experimental 
data for a representative diametric radiation path at x/D = 
0.8 in the flame. Fig. 4-2 shows the representative radiation 
path highlighting the contribution of intensity emitted at 
segment M to the intensity leaving segment N. 

The instantaneous equation of transfer is represented as: 


N [l= 
I, (N) =Liyert, (M4) Tt, (J) (4-6) 
where "N" represents the intensity measurement position and 
"mM" represent the individual radiation segments. This 


formulation requires that the cross correlation between 1T,, and 
I,, be known. In order to treat the effects of turbulence on 
the radiation intensity reaching the detector location "N" 
accurately, simultaneous estimates of emissive power at the 
segment "M" and absorption transmittances at all intervening 
segments "M+1" to "N" would also be required. For the purpose 
of calculating average intensity, estimates of length scales 
and two-point correlation coefficients between I,.(M) and 
(t,.(J), J=M+1,N) are sufficient. However, measurements of 
these quantities are unavailable. To overcome this 
deficiency, the length of the radiation segments must be 
chosen such that it exceeds the local integral length scale. 

In order to study the mean square intensity fluctuations, 
the instantaneous emission intensity for segment M and the 
instantaneous. emission transmittance for segment J are 
transformed following the Reynolds' decomposition more 
commonly used in fluid dynamics: 
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I, (M) -T, (+I. (0) 
t, (J) =ty 3 (0) +04, (J) 


For clarity, the transmittance across the radiation segment 
length is taken as: 


T,e(MN) =[] 7 Tac (D) (4-8) 


After Reynolds' averaging, the mean equation of transfer 


becomes: 
I, (N N) = =Dyeity. (MY) T,, (M, N) (49) 


and that for mean squared intensity fluctuations becomes: 
N 


(K)°(N) = p> (COCs) EMCLErS) ea UCL MC ROLY 
=| 


(4-7) 


+i(M)*(me(M,N))° + 21.)°(M)(4)(M,N) T a(M,N) 


+ 2M) (M)(a.)°(M,N) + 20M) TCM) (RAM, N)) & (MN) 


(4-10) 

Turbulence-radiation interactions in flames originate in 
the non-linearities inherent in the processes of radiative 
transfer (Eqns. 4-4 and 4-5). The non-linearity in the 
equation of transfer is a result of Planck's function (I,,) and 
of monochromatic transmittance (T,) (see® Eqr..975)%. that 
determines the local emission properties. The calculation of 
mean spectral intensity (Eqn. 4-9) is comprised of the mean 
transmittance of local mean spectral emission intensity along 
the radiation path length. Any effects of the fluctuation of 
emission properties on mean spectral intensity would become 
important only if the fluctuations change the mean 
transmittance and mean emission intensity. The calculation of 
the RMS equation of transfer (Eqn. 4-10) includes two 2” 
moment terms, three 3°? moment terms and one 4" moment tern. 
Fluctuating emission properties are seen in every term of this 
equation and should be an important element in the calculation 
of RMS intensity. The equations of the mean and RMS equation 
of transfer (Eqns. 4-9 and 4-10) will be examined to assess 
the validity of the assumption used to overcome the lack of 
spatial correlations for the transmittance and the assumption 
that I',,(M) and T',,(J) are independent. Predictions of mean 
and RMS monochromatic intensity using Eqns. 4-9 and 4-10 and 
local flame emission data are evaluated by comparison with 
monochromatic intensity measurements for segments of diametric 
radiation path shown in Fig. 4-2. 

4.2.3.4 Frequency Content of Emission Properties 

The turbulence radiation interactions originate in 
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Fig. 4-2 Sketch of a radiation path and discretization grid for the equation of radiative 
transfer. 


fluctuating quantities poeeand 1", and, indirectly “in T‘and 
f'\.. These fluctuations have a frequency content determined 
by the local state of turbulence. Fluctuations within a range 
of frequencies are expected to contribute significantly to the 
effects on mean and RMS intensities. For deciding the time 
resolution needed for numerical simulations and experimental 
measurements, knowledge regarding this range of frequencies is 
important. Gaussian filters were used to find the influence 
of frequency content of local transient data on mean and RMS 
intensities for the representative path. 

Data with low frequency content are obtained by using a 
weighted average of neighboring data points. For example, 
intensity at time t is obtained by: 


S=-t+JAt mt 
aT ee (ee em Vc) eet ee) 
where the weight W(J) has a Gaussian form given by: 
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where At = 1/(sampling frequency of the data) and o is the 
standard deviation. The value of o determines the cutoff 
frequency of the filter. 

Fig. 4-3 shows the effect of different o on the frequency 
content of filtered white noise. Normalized power spectral 
densities (PSD) are used to determine the frequency cut off 
level corresponding to each o. MThis cutoff is taken as the 
frequency where the energy content of the PSD begins to drop 
sharply. From Fig.4-3, the frequency cutoffs for o= 3, 10, 
and 15 were found to be 21 Hz, 10 Hz, and 8 Hz respectively. 
Having characterized the filter for different o, its effect on 
turbulent radiation interactions can now be studied. 


4.2.4 Operating Conditions 

The operating conditions are described in Table 2-2. 
Toluene was the fuel used due to its strongly radiating 
nature. The relatively low velocities of the vapor at the 
liquid surface result in complete dominance of the flame by 
buoyancy. The Froude numbers for the present flames at the 
burner surface are higher by over two orders of magnitude 
compared to the acetylene-fired pool flames of Sivathanu et 
al. (1991) and Sivathanu and Gore (1992) allowing a study of 
flame radiation in a different regime. Use of toluene as the 
fuel also separates the effects resulting from peculiar 
chemistry of acetylene from those resulting from strong 
radiative heat loss. The mass burning rates were measured by 
monitoring the output of a load cell mounted under the fuel 
reservoir. The heat release rate is based on complete 
combustion to stable products. The global radiative loss 
fraction is much lower than those of the acetylene/air pool 
flames in spite of comparable soot volume fractions. 
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Fig. 4-3 Power spectral density of filtered white noise 
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4.3 Results and Discussion 

4.3.1 Mean and RMS Properties 

The mean and RMS intensities at 0.9 and 1.0 pum and 
transmittance at 0.632 um are discussed prior to the 
temperatures and soot volume fractions in the flame. 

Fig. 4-4 shows mean and RMS emission intensities at 0.9 
and 1.0 wm as a function of normalized distance from the 
geometric axis of the pool for three positions (x/D= 0.2, 0.4 
and 0.8) relatively close (H,/D = 4.2) to the liquid surface. 
These three locations are in the visually persistent flame 
region. The mean intensities at 1.0 wm are higher than those 
at 0.9 pum consistent with the expected flame temperature 
limits. At x/D =0.2, a narrow flame sheet appears near the 
pool edge (pool edge at r/D=0.5; flame sheet at r/D = 0.42). 
The local emission intensities near the pool center are very 
low for this lowest axial station. The region of high 
intensity is quite narrow at this position. The RMS values 
near the peak intensity are approximately 50% for A = 1.0 um 
and 60 % for A=0.9 wm. These are high RMS values and probably 
originate from the flapping of the flame sheet near the liquid 
surface. As the distance from the liquid surface is increased 
the high intensity region moves closer to the pool center. The 
peak mean intensities at x/D=0.4 are approximately 25% lower 
than those at x/D=0.2. The reduction can be attributed to two 
different reasons, first, the increased RMS values due to 
increased visible flame motion (over 100% for both 0.9 and 1.0 
um) and second the larger time available for heat loss by 
radiation to the surroundings. The region of measurable 
emission intensities broadened by a factor of 2 between these 
two stations. At x/D=0.8, the maximum mean intensity decreases 
to 75% of that at x/D=0.4. For both positions, the RMS 
fluctuations are close to 100% near the center and well over 
100% for the edge locations. The decrease in the maximum mean 
intensity between x/D=0.4 and x/D=0.8 appears to be due to the 
effect of radiative cooling. The width of the measurable mean 
intensity zone now extends across the entire flame width 
possibly signifying a transition from a nearly laminar 
flapping flame sheet at the lower positions to a turbulent 
ensemble of flamelets moving across the entire width of the 
jet. 

: Fig. 4-5 shows the mean and RMS intensities at 0.9 and 
1.0 um for three positions above the persistent flame region. 
Two of the positions are just below the mean visible flame 
height (H,/D = 4.2) while the third is just beyond this 
location. The flame at these farther locations is wider than 
the three positions shown in Fig. 4-4 (note the different 
scale of r/D on the abscissa). For x/D=2.5 and 3.5 positions, 
the peak mean intensities are comparable to those seen at the 
three lower positions. At x/D=4.5, the peak mean intensity has 
decreased substantially, characteristic of the sharp visible 
flame edge exhibited by these flames. The RMS fluctuations are 
much higher at all three positions varying between 100 % and 
250% of the mean intensities. These high RMS values are direct 
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Fig. 4-4 Radial distributions of mean and RMS wideband spectral radiation 
intensity at A = 0.9 and 1.0 wm for locations near the flame base 
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Fig. 4-5 Radial distributions of mean and RMS wideband spectral radiation 
intensity at A = 0.9 and 1.0 um for locations near the flame tip 
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evidence of the effects of turbulent fluctuations on local 
emitted intensities in the highly buoyant flames. 

Fig. 4-6 shows the radial distributions of mean and RMS 
transmittance for three stations (x/D= 0.2, 0.4 and 0.8) in 
the persistent flame region. Consideration of identical 
stations in Figs. 4-4 and 4-6 allows a comparison of the 
similarities and differences in the local absorption and 
emission properties. The mean transmittance at the position 
near the liquid surface is less than unity for positions 
inside of the flame sheet. Visually, a cloud of soot particles 
was observed to move up and down with a low frequency motion 
at this station. These particles do not emit visible radiation 
as seen in Fig. 4-4 and therefore appear to be relatively 
cold. As the region near the flame sheet (r/D = 0.42) is 
approached, the transmittance decreases rather sharply to 
around 0.6 indicating the presence of soot particles. It is 
noted that these particles appear to contribute to the 
emission intensity based on a peak at this location observed 
in Fig. 4-4. Comparison of the transmittance and emission- 
intensity distributions at x/D=0.4 is even more interesting. 
At the first four radial positions near the center at this 
axial station, the mean transmittance based on absorption is 
near 0.5 and relatively constant indicating the presence of a 
large volume fraction of soot. A comparison of these data with 
those of emission intensity (shown in Fig. 4-4) indicates that 
these soot particles are emitting very low intensity light. At 
the radial locations (r/D = 0.35) where the emission 
intensities do increase considerably (change of a factor of 
over 30), the transmittance based on absorption also increases 
to 0.8. This behavior indicates a negative correlation between 
f,, and temperature. The next measurement location (r/D=0.42) 
has an absorption transmittance of 0.95 but still a relatively 
high emission intensity (factor of 10 higher than the central 
region) based on Fig. 4-4. At the last measurement station in 
the persistent flame region (x/D=0.8), the transmittance is 
lowest at the center (similar to that at x/D=0.4) and then 
increases continuously towards the edge. The mean 
transmittance of 0.4 observed at the center-line at this 
position is the lowest measured mean transmittance in this 
flame. Comparison of the radial distribution of mean 
transmittance at this position with that of the visible 
emission intensities (Fig. 4-4) shows that near’ the 
centerline, the lowest intensity occurs with the lowest 
transmittance. The transmittance and the intensity increase 
with radius until the intensity reaches a peak value at r/D= 
0.2. At positions farther than r/D=0.2 from the axis, the mean 
intensity decreases while the mean transmittance continues to 
increase and as expected reaches a value of 1 in the region 
outside the flame. Based on the structure of diffusion flames, 
it is expected that temperature (and emission intensity) would 
reach a maximum at a certain mixture fraction and decrease on 
both fuel rich and fuel lean side. However, whether the 
distributions seen in Figs. 4-4 and 4-6 can be explained by 
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Fig. 4-6 Radial distributions of mean and RMS wideband spectral radiation 
intensity at A = 0.632 um for locations near the flame base 
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the conventional diffusion flame structure (which does not 
account for the strong local radiative heat loss for highly 
sooty regions) is not clear. The RMS fluctuations in 
transmittance are between 0.1 and 0.2 which ranges between 10 
to 50 % of the mean. The reason for these lower values 
compared to the RMS values of emission intensity is discussed 
later. 

Fig. 4-7 shows the radial distribution of transmittance 
for the three axial locations nearer the flame tip. The 
transmittance values are much higher than those seen in Fig. 
4-6 for the positions in the persistent flame region. The mean 
transmittances vary slowly with radial position consistent 
with the emission intensity profiles seen in Fig. 4-5. For 
this region near the flame tip, a comparison of the emission 
intensity and transmittance profiles shows that the emission 
intensity and transmittance change inversely with each other. 
This behavior originates from the higher sensitivity of 
emission intensity to temperature than soot volume fractions 
and is typical of the fuel lean regions in Figs. 4-4 and 4-6 
as well. It is noted that the RMS fluctuations in 
transmittance are approximately 10-25% which are considerably 
lower than the RMS fluctuations in emission intensity (100 - 
300%) seen in Fig. 4-5 for these axial stations. These data 
imply that the emitting material exists in the probe volume 
much more intermittently than the absorbing material. 

The emission intensity measurements summarized above in 
terms of mean and RMS values are interpreted in terms of 
emission temperatures and soot volume fractions based on 
emission in Figs. 4-8 to 4-11. The transmittance data based on 
absorption are interpreted in terms of the soot volume 
fractions in Figs. 4-10 and 4-11. 

Fig. 4-8 shows the radial distribution of mean emission 
temperature at the three axial positions in the persistent 
flame region. In calculating the mean and RMS emission 
temperatures, the readings below the detection limit of the 
PMTs (850 K) have been set to 300 K. This procedure causes a 
downward bias in the mean temperatures and an upward bias in 
the RMS temperatures in the regions of high intermittency in 
the emission signal. 

For radial positions lower than r/D =0.3 at the lowest 
axial station shown in Fig. 4-8, the emission intensity never 
goes above the detection limit of the PMTs. Therefore, this 
region can be classified as relatively cold (at a mean 
temperature of 300 K as per the present technique). The mean 
temperature increases to approximately 1200 K near r/D =0.42 
which is the mean location of the flame sheet. At x/D=0.4, the 
mean temperature near the center is again close to 300 K and 
increases to approximately 1000 K near the flame sheet 
position at r/D=0.35. At this axial station, the radial 
distance traversed by the flame sheet is larger than that at 
r/D=0.2. The RMS temperature fluctuations increase to 40% as 
a result of this motion. Towards the end of the persistent 
flame region at x/D= 0.8, the mean temperature varies between 
900 K and 1100 K for a radial distance of r/D=0.2 before 
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Fig. 4-7 Radial distributions of mean and RMS wideband spectral radiation 
intensity at 4 = 0.632 um for locations near the flame tip 
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Fig. 4-8 Radial distributions of mean and RMS two line emission temperatures 
for locations near the flame base 
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Fig. 4-9 Radial distributions of mean and RMS two line emission temperatures 
for locations near the flame tip 
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Fig. 4-11 Radial distributions of mean and RMS soot volume fractions based on 
emission and absorption data near the flame tip 
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beginning a fast decline in the edge region. The RMS 
fluctuations in temperature vary between 400 K to 500 K. As 
seen in Fig. 4-4, these translate into up to 100% RMS 
variations in the emission intensities. 

The mean and RMS temperature distributions at three axial 
positions nearer to the flame tip are plotted as a function of 
normalized radial position in Fig. 4-7. The temperatures peak 
near the axis for all three positions and decrease 
continuously to the edge. The RMS of temperature fluctuations 
are around 500 K at all three positions. At the highest 
position (x/D =4.5) these happen to be almost 80% of the mean 
temperature level leading to the large variations in 
intensities that were observed in Fig. 4-5. These large 
fluctuations occur due to intermittent behavior of the fire at 
these axial positions. 

Fig. 4-10 shows the radial profiles of mean and RMS soot 
volume fractions based on emission and absorption for the 
three positions closer to the fuel surface. The absorption 
soot volume fractions are a true measure of the volume 
fraction of soot particles in the probe volume while the 
emission soot volume fractions are a measure of volume 
fractions of relatively hot soot particles. If all soot 
particles in the probe volume are at the same temperature, the 
absorption and emission soot volume fractions data should 
collapse within experimental uncertainty limits. Fig. 4-10 
shows that at all positions in the persistent flame region, 
the soot volume fractions based on absorption = are 
substantially larger than those based on emission indicating 
that a large portion of the soot particles are at temperatures 
too low to contribute significantly to the emission intensity. 
At the lowest positions (x/D= 0.2, 0.4), the soot particles 
near the axis seen by the absorption probe (mean of 2 and 11 
ppm) are not seen by the emission probes at all. These 
particles appear as a dark layer near the liquid surface and 
are relatively cold. Presence of large quantities of cold soot 
particles near the surface are expected to affect the heat 
feedback due to blockage. The RMS fluctuations in soot volume 
fractions are typically close to 100 %. Such high fluctuations 
suggest a large level of intermittency in the appearance of 
soot particles associated with soot streaks. 

The radial distributions of mean and RMS soot volume 
fractions based on absorption and those based on emission at 
three positions in the region near the flame tip are shown in 
Fig. 4-11. At all three positions (x/D= 2.5, 3.5 and 4.5) the 
values of soot volume fractions based on absorption are 
substantially higher than those based on emission indicating 
the presence of large amounts of relatively cold soot in these 
regions. The differences between the two measurements 
decrease with axial position probably because the hot soot 
particles that contribute to the emission soot at positions 
near the liquid surface are now at a lower temperature due to 
preferential radiation heat loss. Far away from the flame tip, 
it is expected that the soot particles in the exhaust stream 
reach an isothermal condition with low temperature. Under such 
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conditions the differences between the measurements of soot 
volume fractions based on absorption and emission will 
increase again since f,. Will approach zero and f,, will attain 
a value corresponding to the soot emitted in the exhaust 
stream. The RMS values of soot volume fractions shown in Fig. 
4-11 are comparable to those in Fig. 4-10. This results from 
a combined effect of mixing with ambient air that diffuses the 
soot streaks and the increased turbulence intensity in this 
far-field region. 


4.3.2 Emission Statistics 

4.3.2.1 Flame Properties 

Figs. 4-12 and 4-13 present the character of the local 
emission properties of intensity and temperature at x/D = 0.8. 
The probability density function (PDF) of intensity for three 


radial locations is presented in Fig. 4-12. In the region 
near the pool center (r/D= 0.0 to r/D = 0.14), the PDFs of 
intensity are approximately log-normal. The large dynamic 


range of intensities is made clear by the long tail of the PDF 
curves. The abundance of low intensity points at r/D = 0.28 is 
due to flame intermittency at this position. Fig. 4-13 
presents the PDF of emission temperatures. The detectable 
emission temperatures occur within the range of 1000 to 1700 
K. The large spike at 300 K, especially at r/D = 0.28, is 
indicative of the large amount of low emission intensity 
realizations at this axial position. 

4.3.2.2 Turbulence Radiation Interactions 

Predictions and measurements of mean and RMS intensity (A 
= 1.0 wm) for increasing radiation path length are shown in 
Fig. 4-14. The measured intensity increases rapidly with 
increasing path length due to the increased amount of 
radiating material in the probe volume. At the longer path 
lengths, the intensity levels off due to the effects of self 
absorption and intermittency. 

In order to test the importance of the different 
nonlinear processes, predictions of monochromatic intensity 
for various path lengths were made utilizing: (1) the 
instantaneous measurements of local temperature and soot 
volume fraction and (2) the instantaneous measurements of 
local intensity and soot volume fraction. The results are 
shown in Fig. 4-14. The results indicate that the local 
emission measurements resolved the turbulence scales necessary 
for radiation predictions and that the lack of information 
regarding spatial correlations can be overcome with the 
appropriate choice of radiation segment length. 

Fig. 4-15 presents predictions of intensity which 
utilized the filtered instantaneous intensity and emission 
soot volume fraction measurements. The elimination of 
fluctuations had little effect on the mean intensity levels 
along the path. As turbulence-radiation interactions are 
inherent within the intensity measurements themselves, the 
removal of the high frequency fluctuations would have the 
greatest effect on the non-linear transmittance term. Since 
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Fig. 4-13 PDFs of emission temperature as a function of radial position 
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Fig. 4-14 Predictions and measurements of monochromatic intensity 
(X = 1.0 um) at x/D = 0.8 
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even heavy dampening of the fluctuations had no effect on the 
mean intensity predictions, the non-linearity of the 
transmittance was not significant in radiative transfer. The 
predictions of intensity RMS are affected much more by the 
removal of the high frequency fluctuations. Examination of 
Eqn. 4-10 indicates that fluctuations should play a more 
direct role in the calculation of the RMS of intensity, and 
that is seen by the underprediction of this quantity, even 
when predictions of mean intensity are quite good. 

Fig. 4-16 shows the predictions of intensity using 
filtered instantaneous temperatures and emission soot volume 
fractions. Small reductions in temperature fluctuations (o0= 
3) led to significant underpredictions of both mean and RMS 
intensity. This effect on the mean intensity predictions is 
not seen in Fig. 4-16 and demonstrates the importance of the 
fluctuating quantities when using the nonlinear Planck's 
function for determining the emission intensity term (I,,). 
When the fluctuations of the emission properties are severely 
reduced (o= 15 and 25), the mean and RMS intensity levels drop 


quite dramatically. In fact, the predictions are quite 
similar to those made using mean temperature and emission soot 
properties. 


Fig. 4-17 shows the effect of cutoff frequency on the 
intensity across a 15 mm path length. The mean and RMS 
intensity levels are not affected until approximately a cutoff 
frequency of 11 Hz, which corresponds to the frequency of 
large-scale motion in this flame. 


4.4 CONCLUSIONS 

Use of an existing measurement technique in a liquid 
fired pool flame has resulted in considerable new information 
regarding the mean and RMS distributions of temperature and 
soot volume fractions. These distributions suggest new 
features of strongly radiating flame structure that must be 
considered in the analysis of heat transfer and pollutant 
emission from flames. Predictions of monochromatic emission 
intensity using local emission properties were made. Two 
distinct procedures were used in the treatment of emission 
intensity. The first used local transient emission intensity 
data while the second used local transient temperatures. For 
the unfiltered signals, the results of the two methods are 
identical. However, effects of filtering on the method 
involving temperatures are strong.The specific conclusions 
that can be drawn from the present study are as follows: 


(1) Considerable fluctuations in local emission 
intensities and transmittance measurements based on absorption 
exist in all regions of the present pool flames. 


(Z)) A layer of cold soot particles with large volume 
fractions exists near the liquid surface in the present 
toluene fired pool flame. Based on visual observations these 
particles are probably transported to the cold central region 
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Fig. 4-17 Effect of cutoff frequency on predictions of 
monochromatic intensity (A = 1.0 pm) 
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near the liquid surface by convective velocities and 
thermophoresis. 


(3) The measurements of emission and absorption based 
soot volume fractions in all regions of the flame suggest that 
a considerable portion of the soot particles seen by the laser 
absorption probe are relatively cold and do not contribute to 
emission. The existence of such particles implies a flame 
structure resulting from a balance between radiative heat loss 
and convection and diffusion of energy to regions of high soot 
concentrations. 


(4) Lack of spatial correlations information can be 
overcome with the selection of appropriate radiation length 
scales. 


(5) Turbulence-radiation interactions analysis has 
provided insight into the important physical processes in the 
equation of transfer. A filtering technique involving 


Gaussian distributions along with experimental measurements of 
emission properties was used to verify this analysis. Little 
effect was seen on the mean and RMS intensity levels until all 
frequencies above that of large scale motion for the flame 
were removed. The fluctuations of local emission temperature 
proved to have the most’ significant effect on the 
monochromatic radiation intensity leaving the flame. 


The present measurements indicate that the result of such 
a balance is a relatively low soot layer temperature for the 
heavily sooting toluene flames. The resulting flame structure 
is substantially different from that of weakly radiating 
flames. 
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CHAPTER V 
PREDICTIONS OF RADIATIVE FEEDBACK 


5.1 Introduction 

Studies of radiative heat feedback to the fuel surface in 
liquid pool flames are motivated by the dominant role played 
by this mechanism in the determination of burning rates 
(Burgess and Hertzberg, 1974). This portion of the 
investigation details the local emission/absorption 
measurements in the 30 cm toluene and heptane pool flames. 
The local scalar measurements were used in predictions of 
radiative heat feedback to the fuel surface which will be 
compared with the measurements of this quantity detailed in 
Chapter 3. Many practical fires are heavily sooty and 
measurements of radiative heat feedback under such conditions 
are important. 

Early methods of calculating radiative feedback utilized 
an average flame emissivity, an average shape factor and a 
constant temperature for the entire fire (Burgess and 
Hertzberg, 1974). Multi-ray methods based on the use of 
measurements of local mean absorption coefficients and mean 
temperatures were introduced by Modak (1977). Orloff (1981) 
used Markstein's (1981) scanning radiometer data to obtain 
effective mean temperatures and absorption coefficients as a 
function of height above the pool surface. Variations in 
radiation intensity from different radial positions were 
assumed to originate from different view angles and flame 
shapes. The mean flame shapes were obtained by averaging 
instantaneous photographs similar to the work of Modak (1977). 

Based on measurements of mean and higher moments of 
emission intensity and absorption coefficient for arbitrary 
paths in the pool fires, Markstein (1981) suggested that 
information regarding turbulent fluctuations was necessary for 
obtaining improved models of radiative heat transfer. 

Turbulent fluctuations in temperature and soot 
concentrations are known to cause mean radiative heat fluxes 
to the surroundings to be much larger than estimates based on 
average values (Grosshandler and Joulain, 1986, Gore and 
Faeth, 1986, Sivathanu et al. 1990, Sivathanu 1990). Past 
methods of accounting for turbulent fluctuations have relied 
on stochastic simulations. Some of these studies have used a 
calibrated turbulence model in conjunction with state 
relationships for soot concentrations (Gore and Faeth, 1986), 
while others (Sivathanu et al., 1990, Sivathanu, 1990) have 
avoided these uncertainties by using in-situ absorption and 
emission measurements. 

Sivathanu et al. (1990) and Sivathanu (1990) used two- 
point absorption measurements in conjunction with a soot 
volume fraction-temperature correlation to calculate radiation 
intensity from gas-fired pool flames. The temperature was 
selected between overfire- and underfire- levels based on a 


threshold for emission intensity at one wavelength. The 
results were very sensitive to the selection of the underfire 
temperature. The discrepancies between measurements and 
estimates obtained using low underfire temperatures were 
between 15 and 100%. 

Recently, Sivathanu et al. (1991) have found that a large 
portion of the soot particles measured by absorption probes 
are at relatively low temperatures and do not contribute 
Significantly to the emission of radiation. Therefore, a 
measure of volume fraction of the hot emitting soot particles 
was obtained by utilizing the emission data. A highly negative 
cross-correlation between temperatures and soot volume 
fractions based on emission was observed. A _ bivariate 
stochastic simulation was successfully applied to the 
prediction of spectral radiation intensity using measured 
temperatures (Sivathanu, 1991). 

The stochastic methods of treating turbulence-radiation 
interactions are computer intensive. Therefore, these have 
been applied only to the calculation of directional 
monochromatic radiation intensities leaving the flames along 
few specific directions (Gore and Faeth, 1986, Sivathanu et 
al., 1990, Sivathanu, 1990, Sivathanu et al., 1991). However, 
for engineering applications, estimates of total hemispherical 
intensity and resulting heat flux at points of interest are 
needed. 

Based on this background, the specific objectives of this 
part of the investigation are: (1) to complete local 
measurements of soot volume fractions based on emission 
(representing hot soot particles) and absorption (representing 
all soot particles) and temperatures based on emission in 30cm 
toluene flames; (2) to use new treatment of the equation of 
radiative transfer described in Chapter 4 to account for the 
effects of turbulent fluctuations on the total hemispherical 
heat feedback to the fuel surface; and (3) compare the 
predicted directional and total hemispherical radiative 
feedback with the measurements detailed in Chapter 3. 

5.2 Experimental Methods 

5.2.1 Experimental Apparatus 

The details of the pool burners and fuel supply system 
are described in Chapter 2. All measurements were made in the 
30 cm diameter burner. A constant fuel level of 5 mm below 
the lip of the burner was maintained. 

5.2.2 Instrumentation 

The emission-absorption probe has been described in 
detail in section 4.2.2. The two 6 mm diameter purged tubes 
aligned with 20 mm distance between them defined the 
measurement volume and served as cold background for each 
other. The left side tube collected the emission signal from 
the fire. A beam splitter sent the emission signal to 
photomultipliers equipped with filters with central 
wavelengths of 900 nm and 1000 nm and half-bandwidth of +50 
nm. A second beam splitter guided a laser beam through the 
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measuring volume. The transmitted beam was collected by the 
right hand side purged tube and sent to a laser power meter. 

The signals collected by the calibrated PMTs and the 
laser power meter were amplified and processed through 
antialiasing filters before collection by a _ laboratory 
computer. The signals were sampled at 250 Hz by the computer 
after being filtered at 125 Hz. 

Incident directional radiative heat feedback was measured 
2 mm above the fuel surface using a nitrogen-purged, water 
radiometer. The radiometer and the measurement technique are 
detailed in section 3.2.1. 

5.2.3 Data Analysis 

The methods used to obtain emission temperature, emission 
soot volume fraction and absorption soot volume fraction are 
detailed in section 4.3. If emission intensities were below 
the detection limit of the instrumentation, the emission 
temperature was arbitrarily assigned a value of 300 K and the 
emission soot volume fraction was assigned a value of O ppm. 
The technique used to obtain hemispherical radiative heat 
feedback based on directional intensity measurements is 
described in section 3.3. 

5.2.4 Radiation Predictions 

The radiative heat feedback is obtained by integrating 
directional spectral radiation intensities over the relevant 
wavelength interval and over the hemisphere. As shown in Fig. 
5-1, an imaginary hemisphere on the aperture of the purged 
radiometer is divided into several segments by individual 
rays. Radiation intensity along these rays is calculated 
analogous to the measurements obtained by the narrow angle 
radiometer described in the previous section. 

The calculations involved integration of the equation of 
transfer using discrete radiation segments 1 to N as shown in 
Fig. 5-1 for a representative path. The radiation intensity is 
dominated by soot particles in the present fires. Therefore, 
the mean local spectral intensity emitted by a segment "M" 
along a radiation path is obtained as: 


T,,(M) =-(1-1,,(”) ) Ty, (M) (5-1) 


where I,,(M) is the instantaneous Planck's function defined in 
Eqn. (4-4) and (1-1,,(M)) is the instantaneous emissivity based 
on the transmittance defined by Eqn. (4-5). The overbar 
indicates that average of the product of the Planck's function 
and the emissivity is obtained in order to treat the effects 
of turbulent fluctuations on these non-linear functions of T 
and f,.. The basis for the present calculations is that T(M) 
and f,.(M) are independent of the measurement wavelengths as 
verified by Sivathanu et whet by making measurements at 
several different wavelengths. The K, for all wavelengths 
between 0.5 and 10 um are obtained from Dalzell and Sarofim 
(1969). The fraction of the energy emitted by the segment "M" 
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that is incident on the fuel surface (location "N") can be 
calculated as: 


Tyy(M) -I,,(M TT, Tae (I) (5-2) 


In order to treat the effects of turbulence on the radiative 
heat flux reaching the location "N" accurately, simultaneous 
estimates of emissive power at the segment "M" and absorption 
transmittances at all intervening segments "M+1" to "N" would 
be required. For the purpose of calculating average heat flux, 
estimates of length scales and two-point correlation 
coefficients between I,,(M) and (1,,(J), J=M+1,N) are 
sufficient. However, measurements of these quantities are 
unavailable. It was found earlier that the effects of length 
scales and cross correlations can be treated effectively by 
replacing 1,,(J) by 1,,(J) and then equating the mean intensity 
to the product of the averages of the terms on RHS in Eqn. (5- 
2) (Sivathanu, 1991). As seen from the results, this 
approximation was found to be applicable in the present flames 
indicating that the flame dynamics lead to simultaneous 
occurrence of high I,, (M) and high (1,,(J), J=M+1, N). The 
directional spectral radiation heat flux for a ray can be 
obtained as: 


I, (NM) = ory Tan (ND) (5-3) 


Using these for all rays, the total hemispherical 
radiative heat flux reaching the fuel surface at location "N" 
is calculated as (Siegal and Howell, 1981): 

n/2 f® © x 
q(N) Pi iNet 8 (N) dA)sinOcos6dbdé (5-4) 

The spectral intensities were integrated over 40 segments 
in the wavelength interval 0.5 to 10 um. The resulting 
directional total intensities were integrated over discrete 
rays. The number of rays were varied between 2 and 450 and it 
was found that the total heat flux changes by less than 15% 
between 22 and 450 rays. It is noted that the measurements 
obtained by the gauge correspond to the integral over 
wavelength in Eqn. (5-4). The integration over the hemisphere 
is applied to both calculated and measured intensities to 
obtain the radiative heat feedback. With the above 
formulation, the distributions of mean emissive power and 
emission transmittance are needed for calculating the heat 
feedback including the effects of turbulence. 

5.2.5 Operating Conditions 

The test conditions for local emission/absorption 
measurements of toluene and heptane in the 30 cm burner are 
included in Table 2-2. Radial measurements of local 
emission/absorption properties were made at 10 axial locations 
in the flame in the toluene flame and 8 axial positions in the 
heptane. The axial locations spanned the entire length of the 
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flame. Limited test times (less than 10 minutes) were used to 
control the effects of probe heatup to within 5% of the mean 
emission signal. Measurements along the radial direction 
crossed the entire half width of the flame at each axial 
location. The operating conditions for the radiative feedback 
measurements in toluene are described in section 3.4. 


5.3 Results and Discussion 

5.3.1 Toluene 

Fig. 5-2 shows the radial distributions of the 
measurements of mean and RMS temperatures for 5 axial 
positions. The lowest axial position (x/D = 0.1) is in the 
conical region near the liquid surface. The flame structure at 
this station involves a relatively thin sheet near r/D = 0.55 
where the measured mean temperature reaches approximately 1000 
K. The RMS temperature fluctuations in this region are 400 K 
probably caused by the flapping motion of the flame sheet. At 
all other radial locations (0< r/D < 0.5) at this station, the 
mean emission temperatures are very low as expected for the 
fuel rich locations near the liquid surface. 

At positions farther away from the liquid surface, the 
measurements of mean emission temperatures (see measurements 
for x/D = 0.4 shown in Fig. 5-2 for example) show that the 
structure involves the occasional arrival of high temperature 
material at all radial positions between r/D = 0 and r/D = 
0.25. The mean temperatures are between 1200 K and 1300 K in 
this region with the r/D = 0.2 location showing the highest 
value. The RMS temperatures are near 400 K_ suggesting 
relatively high fluctuation intensities leading to significant 
effects on radiation heat flux. 

Temperature distributions at three axial locations in the 
upper region of the fire are also shown in Fig. 5-2. This 
region is characterized by the highest mean temperature at the 
centerline which decreases due to mixing with ambient air as 
the axial position is increased and very high turbulence 
intensities with the RMS temperature fluctuations approaching 
50% of the mean value. The last position shown in Fig. 5-2 
(x/D = 4.3) has a mean temperature of 450 K at the centerline. 

Fig. 5-3 shows the measurements of emission soot volume 
fractions, f,,, for the same axial positions as in Fig. 5-2. 
Recall that f,, are a good measure of the hot soot particles. 
At the lowest “station, the mean and RMS f,, are negligible due 
to the low temperatures. Near the flame “sheet, the mean f, 
increase to 2 PPM, and the RMS f,, are over 100%. In the region 
beyond the conical zone near the surface, the f,, increase to 
9 PPM and are relatively constant at radial locations between 
r/D=0 and r/D=0.15. As the region of highest mean temperatures 
is approached near the r/D = 0.2 location, the f,, to 
approximately 5 PPM. This is due to oxidation of the soot 
layer in the highest temperature region. Once again the RMS 
fluctuations in the measured f,, are very high suggesting 
significant turbulence-radiation’ “interactions. The 
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Fig. 5-2 Radial distributions of measurements of mean and RMS temperatures 
based on emission at two wavelengths in a 30 cm toluene pool fire 
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Fig. 5-3 Radial distributions of measurements of mean and RMS soot volume 
fractions based on absorption(A = 0.632 um) in a 30 cm toluene pool fire 
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measurements in the upper region show that the f,, decrease 
continuously as the temperature decreases. Finally, at the x/D 
=4.3 location, the mean f,, approach 0.2 PPM with over 400% RMS 
fluctuations. 

Radial distributions of soot volume fractions based on 
absorption, f,,, are shown in Fig. 5-4 for the same 5 axial 
locations. At the lowest position the mean f,, reach their peak 
value (10 PPM) at the center and remain relatively high in the 
core region up to r/D=0.2 and then slowly decrease near the 
flame sheet to approximately 6 PPM. Visual observations 
indicate that soot particles are carried to the central cold 
region near the fuel surface by the complex hydrodynamics of 
the entrainment and fire induced flow. The existence of cold 
soot particles near the liquid surface probably leads to 
substantial blockage of radiative heat feedback. It is noted 
that even at the location of highest mean temperature, the 
measurements of f,, are approximately 3 times those based on 
emission suggesting the presence of relatively cold soot near 
the flame sheet. The RMS fie are very high even in the region 
near the fuel surface similar to the observations of Bouhafid 
et al) afi989)- 

The mean f,, are higher than f,, at almost all the axial 
stations. At x/D=2. 3, the mean f,, and f,, are close suggesting 
local mean homogeneity of the material in the probe volume at 
this axial station. As this material cools by mixing with 
Surrounding air and by radiative heat loss, f,, exceed f,, 
substantially. Finally, the soot particles observed by the 
absorption probe at x/D=4.3 represent the thick dark smoke 
that is visually observed in the region above the fire. 

The measurements and calculations of directional total 
radiative heat flux are plotted as a function of radial 
position on the pool surface for three representative rays in 
Fig. 5-5. The ray marked 8 = 90°, #@ = 0° is in the direction 
of gravity and shows the highest measured heat flux. The ray 
marked 86 = 110°, # = 0° is directed at 20° to the direction of 
gravity and starts away from the fire axis. The ray marked 6 
= 130°, @ = 0° is directed at 40° to the direction of gravity 
and starts away from the fire axis. Due to the shorter path 
length in the fire, this ray shows the lowest heat flux. The 
calculated directional heat fluxes are in reasonable agreement 
with the measurements considering the uncertainties in soot 
refractive indices, length scales and spatial cross 
correlations in soot concentrations and temperature. 
Discrepancies larger than the uncertainties occur for the ray 
in the direction of gravity at r/D=0.2. However, the overall 
performance of the radiation analysis is encouraging. 

Fig. 5-6 shows measurements and predictions of total 
hemispherical radiation heat flux to the surface of the 
toluene pool obtained by integrating directional fluxes 
Similar to those shown in Fig. 5-5. The data were obtained 
from the integration of 22 directional fluxes with appropriate 
weights. Calculations based on 22 to 450 rays were performed 
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Fig. 5-5 Measurements and calculations of hemispherical radiative intensity incident 
the liquid surface as a function of radial location and view angle 
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Fig. 5-6 Measurements and calculations of hemispherical radiative heat flux incident 
the surface of a 30cm toluene pool fire as a function of radial location 
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and showed 15% variation. The results reported in Fig. 5-6 are 
based on 450 rays. The agreement between measurements and 
calculations of heat feedback is good supporting the present 
method. Based on an energy balance at the liquid surface, the 
radiative heat flux of approximately 30 kW/m? is sufficient to 
maintain the observed burning rate. Thus the fire is dominated 
completely by radiative heat feedback. 


5.3.2 Heptane 

Measurements of the absorption/emission properties 
in heptane show a different flame structure than that seen in 
toluene. Temperatures in the lower region of the flame (x/D 
= 0.2) as seen in Fig. 5-7 show that mean temperatures are 
much higher than seen in the toluene flame at a similar axial 
position. Even at higher axial locations (e.g. x/D = 2.3) the 
mean temperatures are generally greater than found in the 
toluene flame. 

Fig. 5-8 presents the radial distribution of emission 
soot volume fractions at three different axial locations. As 
seen in the previous flames, the RMS values of this quantity 
are quite large, sometimes double that of the mean soot volume 
fraction. The largest amount of emission soot occurs at a 
much higher location in the heptane flame (x/D = 1.5) as 
compared to the toluene flame (x/D = 0.4). 

Fig. 5-9 shows measurements and predictions of total 
hemispherical radiation heat flux to the surface of the 
toluene pool obtained by integrating directional fluxes. The 
calculations are based only on soot emission and neglect the 
contribution of gas band radiation, which should be 
considerable in the heptane flame. Calculations of radiative 
heat flux to region near the burner edge (r/D > 0.2) were 
quite low due to the coarseness of the grid of 
emission/absorption measurements, especially near the fuel 
surface. 

Fig. 5-10 shows the contribution of radiant heat flux to 
the pool center from various axial distances from the fuel 
surface based on predictions. The results indicate that only 
radiation from the region near the fuel surface (x/D < 1.0) 
contributed to the feedback. Radiation from regions above x/D 
= 1.0 does not reach the fuel surface due to absorption. The 
contribution of radiation heat flux from different axial 
locations in heptane is also predicted. Measurements of the 
emission and absorption properties were made in the manner as 
described in Section 5-3. The predictions of heptane feedback 
are only based on soot emission, though gas emission is 
thought to be important. Contributions to the radiation 
incident on the fuel surface are seen from most axial heights 
in heptane though the largest contribution does originate from 
near the fuel surface. The self absorption seen in toluene 
is not as important in the heptane flame, due to lower amount 
of soot present. 
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Fig. 5-7 Radial distributions of measurements of mean and RMS 
temperature based on emission at two wavelengths in a 30 


cm heptane pool fire 
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Fig. 5-8 Radial distributions of measurements of mean and RMS 
soot volume fractions based on emission at two 
wavelengths in a 30 cm heptane pool fire 
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Fig. 5-9 Measurements and calculations of hemispherical radiative 
heat flux incidient on the surface of a 30 cm heptane pool 
fire as a function of radial location 
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5.4 Conclusions 

Detailed emission and absorption measurements as well as 
directional and hemispherical total heat feedback measurements 
in 30 cm toluene fires have shown that: 


(1) Large fluctuation intensities in temperatures and 
soot volume fractions based on absorption and emission exist 
in all regions of the fire. This is similar to the structure 
seen in the toluene flame for the 7.1 cm burner (Chapter 4) 
and that observed in ealier studies (Bouhafid et al., 1989; 
Fischer et al., 1987). 


(2) A transient layer of cold soot particles exists near 
the liquid surface. The role of this layer in blocking the 
radiative heat feedback can be treated within the present 
radiation analysis. 


(3) A significant fraction of soot particles at most 
locations in the flames are at relatively low temperatures and 
do not contribute to the radiative feedback. 


(4) The radiation analysis based on _- transient 
measurements of temperature and soot volume fractions can be 
applied successfully in conjunction with ae simplified 
treatment of turbulence radiation interactions to estimate 
directional total and hemispherical total heat feedback to the 
fuel surface. However, predictive models for obtaining the 
necessary transient distributions are sparse. 
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CHAPTER VI 
SUMMARY AND CONCLUSIONS 


6.1 Summary 

The overall objective of the present study was to the 
improve the understanding of radiation from liquid fueled pool 
flames. Particular emphasis was placed on investigating the 
mechanism of radiative heat feedback in pool flames. The 
effect of fuel type and scale on burning rates and flame 
radiation was studied for a variety of alcohol, parafin, and 
aromatic fuels. A multi-location method was used to measure 
total radiative output from these flames. Measurements were 
made of the radiative energy incident on the fuel surface. 
Local emission and absorption properties were determined for 


all regions of luminous pool fires. Predications of 
directional and hemispherical radiative heat feedback were 
made using the local emission/absorption properties. The 


effects of turbulence-radiation interactions on radiation 
predictions of radiative heat feedback and monochromatic 
intensity were investigated. 


6.2 Conclusions 

In addition to the conclusions presented at the end of 
each section, the significant observations and contributions 
of the present study are as follows: 


al The total radiative heat loss fractions for a 
particular fuel were found to be constant for the scales 
investigated. Families of fuels were also found to have 
comparable radiative heat loss fractions. 


2. Measurements of radiative heat feedback showed that 
radiation is a significant heat transfer mechanism in sooting 
flames but also in non-luminous flames. Radiation contributed 
almost the entire heat transfer to the fuel surface in toluene 
and heptane flames. Radiative heat feedback was significant 
in the center of methanol flames. 


& is The emission/absorption properties of pool flames 
were found to have a large dynamic range. Significant 
fluctuations of scalars were found in all regions of the 
flames. Large amounts of non-emitting soot were found in all 
regions of a strongly radiating pool flame. This finding 
explains the similarity in radiative loss fractions for fuels 
of different sooting tendencies. 


4. The importance of turbulence-radiation interactions in 
radiative heat transfer was shown. Fluctuations in scalar 
properties were shown to be important for predictions of 
radiative heat feedback and monochromatic intensities. 


6.3 Recommendations for Further Research 
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The present study concentrated on the local emission 
properties of soot dominated, strongly radiating flames. [In 
optically thin flames, radiation from combustion gases must be 
considered. Simultaneous measurements of radiating gas 
concentrations in pool flames must be made. Planar 
measurements of emission properties, as opposed to the single 
point measurements made in this study, would give additional 
information about turbulence-radiation interactions in pool 
flames. 
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HEPTANE 7.1cm TOTAL RADIATIVE OUTPUT CALIBRATION CONST = 0.208 W/cm2-mv 
hep2vert.dat hep2rad.dat 
axial radial 


-36820E+00 6.35 


.27100E+00 8.89 

-20090E+00 11.43 
-79360E-01 10.16 .16200E+00 13.97 
.82290E-01 15.24 -14600E+00 16.51 
-79000E-01 17.78 -110S50E+00 19.05 


.69750E-01 30.48 -46230E-O01 31.75 
-65350E-01 _ 33.02 .34980E-01 34.29 


SIDE FLUX= 249.0908W BOTTOM FLUX=189.1772W 


TOTAL FLUX= 4.382681E-01 kW 
HEPTANE 30CM TOTAL RADIATIVE OUTPUT 


HAX12.DAT HRA12.DAT 
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HEPTANE 100cm TOTAL RADIATIVE OUTPUT 
hep#6xrc.dat hep#6xrr.dat 
AXIAL RADIAL 


Q, W/cm2 POS (cm) Q, W/cm2 POS (cm) 


0.36700E+00 0.00 0.33340E+01 0.00 


0.42638E+00 50.00 0.33340E+01 60.00 


0.44623E+00 100.00 0.15395E+01 122.00 


0.46634E+00 150.00 0.31405E+00 250.00 


0.36206E+00 290.00 0.50996E-01 345.00 


SIDE FLUX= 470250W BOTTOM FLUX= 164269.0W 


TOTAL FLUX= 634.519kW 
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HEPTANE 30CM TOTAL RADIATIVE OUTPUT 
HAX12.DAT HRA12.DAT 
axial radial 


Q, (W/cm2) POS (cm) Q, (W/cm2) POS (cm) 


.38039E+00 .00 .91546E+00 00 


SIDE FLUX= 35479.58W BOTTOM FLUX= 3488.017W 


TOTAL FLUX= 38.9676KW 
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HEXANOL 7.1cm TOTAL RADIATIVE OUTPUT CALIB. CONST = .208 W/cm2-mv 


HEXALCVT.RAD HEXALC1B.RAD 


POS (cm) SIGNAL (v) 


.35450E-01 .00 .28490E +00 .00 
.35340E-O01 2.54 .28490E+00 6.35 
.38220E-01 5.08 -18860E+00 8.89 
.37120E-01 7.62 -12430E+00 11.43 
.38080E-01 10.16 -87020E-01 13.97 
.36570E-O1 12.70 -67150E-01 16.51 
36110E-O1 15.24 -49440E-01 19.05 
.37820E-O01 17.78 .36000E-01 21.59 


.36920E-01 20.32 -29190E-01 24.13 


-32100E-O1 22.86 .22880E-01 26.67 
.29810E-O1 27.94 -19380E-01 29.21 
-26800E-01 33.02 -15690E-O1 31.75 
-20470E-01 38.10 -11410E-01 36.82 
.15620E-O1 43.18 -508S0E-02 41.91 


-91480E-02 53.34 


.00000E+00 73.61 


SIDE FLUX= = 214.788800W BOTTOM FLUX= 84.783150W 
TOTAL FLUX= 2.995719E-OIkW 


MASS FLOW RATE= 2.32g/min 
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METHANOL 7.1cm TOTAL RADIATIVE OUTPUT 


METFLUXV.AVG METFLXRB.AVG 
axial radial 


SIGNAL (v) POS (cm) SIGNAL (v) POS (cm) 


-40920E-01 .00 .19296E +00 .00 
44989EO1 2.54 -19296E+00 6.35 
-4737SE-O01 5.08 -11949E+00 8.89 
SOSOSE-O1 7.62 -80454E-01 11.43 
-50366E-O1 10.16 59800E-01 13.97 
-S160S5E-01 12.70 -44879E-01 16.51 
.52936E-O1 15.24 .36288E-01 19.05 
.52350E-01 17.78 -28638E-O01 21.59 
.45364E-01 20.32 .24460E-01 24.13 
.44866E-O1 22.86 .20701E-01 26.67 


-4S9SIE-O1 25.40 -I817S5E-O1 29.21 
.37820E-01 27.94 
.35079E-01 30.48 
.29864E-01 33.02 
.28021E-01 35.56 
.28443E-01 38.10 
.25831E-O1 40.64 
.23932E-01 43.18 
-20620E-O01 45.72 
.17984E-O1 48.26 
-16657E-01 50.80 


.00000E+00 67.55 


SIDE FLUX= 189.4828W BOTTOM FLUX= 48.27W 
TOTAL FLUX= 2.377575EO1kW 


RADIATION FRACTION= 2.0E-01 MASS FLOW RATE= 3.54g/min 
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METHANOL 30cm TOTAL RADIATIVE OUTPUT 


m2ax.dat m2ra.dat 
axial radial 


Q, (W/cm?) POS (cm) Q, (W/cm2) POS (cm) 


.56400E-01 .00 -26871E+00 .00 
.60967E-01 5.00 -26871E+00 33.65 
.62800E-01 10.00 .17226E+00 38.65 
-65772E-01 15.00 -12019E+00 43.65 
-63060E-01 20.00 .87388E-01 48.65 
-65616E-O1 25.00 -63150E-01 52.65 
-62014E-01 30.00 -47S510E-O01 57.65 
-61195E-01 35.00 -36852E-01 62.65 
-62186E-O1 40.00 .28524E-01 67.65 
-ST0SIE-O1 45.00 .24192E-01 72.65 
.53792E-01 50.00 

.48065E-01 60.00 

.39170E-01 70.00 

.33486E-01 80.00 

.26143E-01 90.00 

.22318E-01 100.00 

-17692E-01 110.00 

-15746E-01 120.00 


-00000E+00 162.17 


SIDE FLUX= 2535.563W BOTTOM FLUX=876.90W 
TOTAL FLUX= 3.412463kW 


RADIATION FRACTION= 0.17 MASS FLOW RATE= 0.9 g/s 
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METHANOL 100cm TOTAL RADIATIVE OUTPUT 
met#8xrc.dat met#8xrr.dat 
’ AXIAL RADIAL 

Q, W/cm2 POS (cm) Q, W/cm2 POS (cm) 


0.22973E-01 0.00 0.41005SE+00 0.00 


0.35090E-01 50.00 0.41005E+00 60.00 


0.32085E-01 100.00 0.46911E-01 125.00 
0.31711E-O1 150.00 0.15743E-01 192.00 
0.22722E-01 222.00 0.68980E-02 230.00 
0.16869E-O1 292.00 0.11580E-02 275.00 


0.00000E+00 446.00 0.65500E-03 330.00 


SIDE FLUX= 19418.5W BOTTOM FLUX= 8024.07W 


TOTAL FLUX= 27.4425kW 
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TOLUENE 7.1cm TOTAL RADIATIVE OUTPUT CALIB. CONST = 0.208 W/cm2-mv 


tol2vert.ave tol2rad.ave 
axial radial 


SIGNAL (v) POS (cm) SIGNAL (v) POS (cm) 


18150E+00  .00 .14390E+01  .00 
19440E+00 2.54 14390E+01 6.35 
.19020E+00 5.08 .84520E+00 8.89 
19930E+00 7.62 53660E+00 11.43 
20990E+00 10.16 39880E+00 13.97 
20550E+00 12.70 32060E+00 16.51 
21190E+00 15.24 25450E+00 19.05 
21620E+00 17.78 20440E+00 21.59 
20980E+00 20.32 14940E+00 26.67 
20630E+00 22.86 10070E+00 31.75 
.19900E+00 25.40 72190E-01 38.10 
18830E+00 27.94 66490E-01 39.37 
16750E+00 33.02 
15180E+00 38.10 

" 13440E+00 43.18 
.11500E+00 48.26 
97830E-01 53.34 
84510E-01 58.42 
63910E-01 66.04 


-00000E+00 88.31 


SIDE FLUX= 1415.985W BOTTOM FLUX= 448.3577W 
TOTAL FLUX= 1.864343kW 


RADIATION FRACTION= 0.312 MASS FLOW RATE= 8.73 g/min 
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TOLUENE 100CM TOTAL RADIATIVE OUTPUT 
tol#4xrc.dat tol#4xrr.dat 
AXIAL RADIAL 


Q, W/cm2 POS (cm) Q, W/cm2 POS (cm) 


0.33871E+00 0.00 0.41444E+01 0.00 
0.39228E+00 50.00 0.41444E+01 60.00 
0.40000E+00 100.00 0.14938E+01 122.00 
0.40878E+00 150.00 0.13703E+00 250.00 
0.33666E+00 222.00 0.44808E-01 285.00 
0.19149E+00 350.00 0.19740E-01 345.00 


0.11855E+00 400.00 0.84190E-02 442.00 


0.00000E+00 502.00 


SIDE FLUX= 366777W BOTTOM FLUX= 154290 W 


TOTAL FLUX= 521.067 kW 


ay, 


APPENDIX C 


DIRECTIONAL RADIATIVE FEEDBACK 
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APPENDIX D 


EMISSION/ABSORPTION PROPERTIES 
TOLUENE 30 cm 
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APPENDIZ E 
SPECTRAL EMISSION AND FUEL VAPOR TRANSMITTANCE 
FUEL SURFACE REFLECTION 
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FLAME EMISSION 


TOLUENE HEPTANE 


A I A I 
(nm) (W/m?-sr-um) (nm) (W/m?-sr-pum) 
850.0 144.0 1200.0 697.8 
900.0 206.3 1300.0 766.4 


950.0 237.2 1400.0 847.0 
1000.0 329.7 1500.0 962.7 


1600.0 2664.3 2000.0 898.1 
2098.8 2100.0 517.7 
1800.0 2897.3 2200.0 609.4 


2300.0 2887.0 2700.0 1078.1 
2400.0 2990.0 2800.0 1273.5 


3200.0 1480.0 3500.0 892.7 
3400.0 1634.0 3600.0 505.2 


4400.0 3186.0 4200.0 1128.8 


1700.0 
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4600.0 1789.7 4300.0 781.5 


| 4800.0 1739.0 4400.0 3950.1 


4500.0 2560.6 


4600.0 1601.6 


4700.0 409.9 
4800.0 248.5 


FUEL VAPOR TRANSMITTANCE 


TOLUENE HEPTANE 


a TRANSMITTANCE a TRANSMITTANCE 
(nm) (nm) 


1000.000 0.940 1200.000 0.955 
1100.000 0.947 1300.000 1.000 
1200.000 0.965 1500.000 0.997 
1300.000 0.957 1600.000 1.000 
1400.000 0.959 1700.000 1.000 
1500.000 0.959 1800.000 0.998 
1600.000 0.965 1900.000 1.000 


1700.000 0.950 2000.000 0.995 
2100.000 0.997, 


0.972 


1800.000 


1900.000 0.980 2200.000 0.998 


2000.000 0.969 2300.000 1.000 


2100.000 0.968 2400.000 0.957 


2200.000 0.946 2500.000 0.995 
2300.000 0.945 2600.000 0.889 
2400.000 0.941 0.965 
2500.000 0.915 2800.000 0.931 
2700.000 0.963 3000.000 0.950 
2800.000 0.960 3100.000 0.940 
2900.000 0.967 3200.000 0.918 
0.996 


2700.000 


3300.000 


3000.000 0.955 


132 


4000.000 0.999 3800.000 0.936 
4200.000 0.989 3900.000 0.956 


_ 2 sine reel re 


FUEL SURFACE REFLECTION 


HEPTANE 


TOLUENE METHANOL 


WATER 


ANGLE REFLEC. ANGLE REFLEC ANGLE REFLEC. ANGLE REFLEC. 


6.49 0.580 6.49 0.660 5.37 0.444 6-60970.925 
11.00 0.370 12.10 0.360 13.20 0.327 14.50 0.263 
18.00 0.130 26.40 0.126 2225900800017 22.30 0.120 


35.50 0.074 35.50 0.075 54.78 0.026 54.78 0.036 
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